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LXXVIIT. Simultaneous Diffusion and Reversible Chemical Reaction 


By J. Crank 
Courtaulds Limited, Maidenhead, Berks.* 


[Received March 21, 1952] 


ABSTRACT 


The problem considered is one of diffusion in which some of the diffusing 
substance is immobilized by chemical reaction as diffusion proceeds. The 
reaction is considered to be first order and reversible, the forward reaction 
proceeding at a rate proportional to the concentration of solute free to 
diffuse, and the backward reaction at a rate proportional to the concen- 
tration of immobilized solute. Mathematical solutions are derived for 
the diffusion of a limited amount of solute into a plane sheet, a cylinder, 
and a sphere respectively. The solutions comprise two infinite series and 
the physical significance of the various terms is discussed briefly. Solu- 
tions involve three independent parameters ; these are the ratio of the 
volumes occupied by solution and sheet, cylinder or sphere respectively, 
the partition factor between immobilized and free solute in the final 
equilibrium state, and a modified rate constant expressing the relative 
rates of diffusion and reaction. Calculated values of M,/M,, where M, 
is the total amount of solute present in the sheet, cylinder or sphere at 
time ¢t and ,, the corresponding quantity after infinite time, are tabulated 
for the case of an infinite amount of solute and for a range of values of the 
other two parameters. Characteristic features of diffusion-with-reaction 
are shown graphically. An important general conclusion is that if the 
diffusion is more than a thousand times faster than the reaction, expressed 
in terms of the respective half-times of the two processes, the behaviour 
of the joint diffusion-reaction process is almost the same as if diffusion 
were infinitely rapid. If on the other hand the half-times for diffusion and 
reaction are comparable the behaviour approximates to that for an 
infinitely rapid reaction. It is shown that the mathematical solutions 
derived for absorption also describe the course of desorption under 
comparable conditions. Solutions for an irreversible reaction are deduced 
as special cases of those for a reversible reaction. 


$1. INTRODUCTION 


THE problem to be discussed is that of the absorption of one substance by 
another through which it can diffuse and with which it can also react 
chemically. This can be regarded as a problem in diffusion, in which some 
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of the diffusing substance becomes immobilized as diffusion proceeds, or as 
a problem in chemical kinetics in which the rate of reaction depends on the 
rate of supply of one of the reactants by diffusion. There are numerous 
practical examples of processes involving simultaneous diffusion and 
chemical reaction of one sort or another. Thus diffusion may take place 
within the pores of a solid body which can absorb some of the diffusing 
substance; or we may have diffusion occurring through a gel and an 
immobile product resulting from the attraction of the diffusing molecules 
to fixed sites within the medium. This occurs, for example, when sub- 
stantive dye molecules diffuse into water-swollen cellulose. Similar 
processes may be involved in the swelling of a high polymer substance by a 
solvent. In other cases, the diffusing substance may be immobilized by 
the formation of anew compound between itselfandthemedium. Examples 
involving diffusion into living cells and micro-organisms can be cited from 
biology and biochemistry. One example from physiology which has been 
considered in detail recently by Nicolson and Roughton (1951) is the 
uptake of oxygen by red blood corpuscles. Chemical reactions in high 
polymer substances are often considerably dependent on the mobility 
of the reactants as well as on the kinetics of the reaction itself. 

It this paper we confine attention to cases in which the process by which 
the diffusing substance is immobilized is first order and reversible. This 
appears to be the most general case for which formal mathematical 
solutions can be obtained at the present time. Although of course in 
practice not all processes will be strictly first order, nevertheless it is 
thought that the number for which it will be a useful approximation is large 
enough to justify extensive mathematical treatment. Subject to this 
condition the treatment is quite general, irrespective of the precise nature 
and mechanism of the diffusion or of the formation of immobilized product. 
The aim of this paper is not to examine any one practical example in 
detail. It is rather to establish quantitatively the characteristic features 
of the diffusion-with-reaction system, and to provide theoretical curves 
which the experimentalist can use in attempting to assess his practical 
Measurements. 

The behaviour to be expected when the reaction is reversible clearly 
depends on the relative rates of diffusion and reaction. Several authors 
(Hill 1929 ; Wilson 1948 ; Crank 1948) have discussed problems in which 
the product is formed so rapidly compared with the diffusion that it is 
always effectively in equilibrium with the diffusing substance. At the 
other extreme is the case of diffusion being so rapid compared with the 
reaction that the concentrations of diffusing substance and product are 
effectively uniform throughout the medium and the behaviour is controlled 
solely by the reversible reaction. So far as is known only one mathematical 
solution has been obtained (Wilson 1948) for the more general case in which 
the rates of diffusion and reaction are comparable. In attempting to 
evaluate Wilson’s series solution numerically some interesting features 
concerning its physical significance, and also some unexpected difficulties 
in evaluation emerged, and this led to the problem being taken up anew. 
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§2. STATEMENT OF THE PROBLEM 

The general problem considered in this paper can conveniently be stated 
in terms of a solute diffusing from a solution into a plane sheet of material. 
The modifications necessary for corresponding, alternative problems, such 
as those of a sphere or cylinder suspended in a vapour, will be obvious. 

Suppose an infinite sheet of uniform material of thickness 2a is placed in 
a solution and that the solute is allowed to diffuse into the sheet. As. 
diffusion proceeds, a first order, reversible reaction occurs and a product, 
which is non-diffusing, is formed. The sheet occupies the space 
—a <x <a, and there is a restricted amount of solution which occupies the. 
space —l—a<x<—a, a<u<l+a. The concentration of solute in the- 
solution is always uniform and is initially Cy, while initially the sheet is. 
free from solute. Let c be the concentration of solute free to diffuse 
within the sheet and s that of the immobilized solute, each being expressed. 
as amount per unit volume of sheet. 

The diffusion is governed by the equation 


dc dc ods 
Ape Gate Fe (2.1) 
and we consider the simultaneous reaction to be of the type 
GS Oba NCS, Oar a oe aps be cs 3 Bos. (222) 


Here DP is the diffusion coefficient and A and p are the rate constants of the 
forward and backward reactions respectively. Thus the immobilized 
solute is formed at a rate proportional to the concentration of solute free 
to diffuse, and disappears at a rate proportional to its own concentration. 
We require solutions of (2.1) and (2.2) with the initial condition‘ 

SG) eee One tases pve os, cw (200) 
and with a boundary condition expressing the fact that the rate at which. 
solute leaves the solution is equal to that at which it enters the sheet over 
the surfaces x= +a. This condition is 

loci0t—=—- Docjox, w= 0,6 >0.° — 2 2. (BRED 
We here assume that the concentration of solute free to diffuse just within. 
the surface of the sheet is the same as that in the solution. This may not: 
be so and there may be a distribution factor K, which is not unity, such 
that the concentration just within the sheet is K times that in the solution. 
This can clearly be allowed for by using a modified length of solution, 
1/K, in place of J in (2.4) and elsewhere. Mathematical solutions follow 
for these equations and for corresponding equations for the cylinder and 
sphere. 

§3. MarTHEMATICAL SOLUTIONS 
(a) Plane Sheet 

Solutions of the equations of § 2 can be obtained by the method of 

Laplace transforms (Carslaw and Jaeger 1941). Writing ¢ and § for the 
Laplace transforms of c and s respectively, so that 


¢= if c exp (— pt) dt, 5 | s exp (—pft) dt, a ey 
0 0 


epee 
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we have the following equations for ¢ and s: 


pé=—ps+Dorc/oa?, . . . . » «, (3-2) 
pi AE— ps ,) oe Se ee ae 
—IC,+plé=—Doéjdx, =a. 2 1 2 ee (Sa) 


On eliminating § from (3.2) and (3.3) and replacing the partial derivative 
by an ordinary derivative since ¢ does not appear, we find 


de z Pp ptate 
zs 2¢—0, P®=—=———_, . . . . . (3.5 
dx2 +k Cc 0, k D pte ’ ( ) 
of which the solution that gives é an even function of x is, 
é==F (p) cos kz... 0 oe 


The function F(p) is determined by the boundary condition (3.4) and it 
follows immediately that 
IC) cos kx 


ira pl cos ka—kD sin ka’ 


(3.7) 


The application of the Inversion Theorem (Carslaw and Jaeger 1941) is — 
straightforward and after some reduction gives 


IC as Co exp (Prt) cos k,,x 
‘ae I+-(R-+1)a + 4 Au a De p,’la | cos k,a’ 
MM Giermel Lat aks aah} 
F (3.8) 
where the p,,’s are the non-zero roots of 
l 
oa =k, tank,a, k,?2=— ee : (3.9) 


and R=Ajp is the partition factor between immobilized and free solute. 
The expression for s differs from (3.8) only by having an. extra factor 
A/(p, +) multiplying the nth term, including the term p,=0. Writing M, 
for the total amount of solute, both free to diffuse and immobilized, in 
half the sheet at time ¢, and ,, for the corresponding quantity in the final 
equilibrium state attained theoretically after infinite time, we have 


M (oe) 
M., ty ty4 Au a, Pn Pyila 
(p,teef 120 Dk * BDF 
where the p,,’s are given by (3.9) and where 
a=I/(R+1)a, M,=IC,/(1+a). oo tua a eee 


(6) Cylinder 
The case of the cylinder has previously been considered by Wilson (1948), 
using a slightly different method. The final result for a cylinder of radius, 
a, in a solution occupying a region of cross-sectional area, A, is 
My (1a) exp (Ppt) 
MM i a Au ma? AD AANS 
Ih 4l+ ot + oe 
@,+epef Va? DES 


(3.12) 
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where 
We Aran ives) A Me AG | (l-po)) a a 2 ar (3.13) 
and the p,’s and k,’s are given by 


Ap, "se. J (kn) Et _ PnPataAth 


27Da~ "J (ka)? OE SO re 
The expression for the concentration of free solute is 
aa ACy > c EXP (Prt) J (knit) 
oo A-+(R--1)ma? ers Co eee ma* = Ap,? ) Jo(kya) gn) 
oat ae tebe 


The expression for s differs from (3.15) only by having an extra factor 
A\/(P,+) multiplying the nth term, including the term p,—0. 
(c) Sphere 

On introducing new variables rc and rs the equations for the spherical 


case take essentially the same form as those for the plane case, and the 
same method of solution leads to 


M. 1y in Au 27a Pepe ye 
seh ele (posta EY. =e ic saypees 
where the p,’s and k,,’s are given by 
er sy, 2 PaPat Ate 
TnDa 2” l= 0 Cot ka, 3k, — Dae. (3.17) 
and 
3V VG, 
~ 47ra3(R+1)’ HES Ieee: CRE) 
The expression for c is 
Ise 3VC, 
rag 30+ na R+1) 
ag x Co exp (Dut ) a sim ken (3.19) 


ma f paces eo Hae Pn Pine Pane Vp, )rsink,@’ 
Bat eee — 3Dk,2 * BnaDh,2 


and that for s differs from ae only by having an extra factor A/(p,+) 
multiplying the nth term, including the term p,=0. 


§4. PHysicaL SIGNIFICANCE OF THE MATHEMATICAL SOLUTIONS 
When mathematical solutions are as complicated in form as those 
in § 3, their physical significance is not immediately obvious. Consider 
eqns. (3.9) and put 


c= ka, y=pa?/D, €=a(A+p)/D,n=a'p/D. . . . (4.1) 

Then 
iG ane tan) Gh ee oe ee  E.8) 
spe a nip a1) 


yn’ 
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are the equations to be solved for the roots. Graphs of pede a 
are sketched in fig. 1 to show the general location of the roots. e grap 
of (4.2) is the same for all € and 7, and from the figure it is easy to ps: 
qualitatively how the roots vary with € and n. When pees a 
transformations to those of (4.1) are applied to eqns. (3.14) for the 
cylinder and (3.17) for the sphere, the resulting equations are of the ae 
form as (4.2) and (4.3), so that fig. 1 can be taken as showing qualitative y 
the location of the roots for all three cases. There is a root for which 
k,a is imaginary given by 

lylja=/(—2) tanh +/(—2) . oe Oa) ee (Soca) 


calculated from (42) y 
—- eee calculated from (4:3) 
1O successive real roots 
@ ‘oot for which 2,2 is imaginary 


ord. | y abs. | a 


Location of roots of (4.2) and (4.3). 


‘The general expression for M,/M,, therefore comprises a unit term from 
Pn—9, a term for which k,a@ is imaginary, and two infinite series of terms 
corresponding to the intersections of the two branches of (4.3) with 
successive branches of (4.2). 

The relative importance of the various terms depends on the parameter 
n. It is interesting to see quantitatively what happens to the general 
solution for the extreme values of this parameter which correspond to very 
fast and very slow reactions. The roots of (4.3) are given by 


2y=—(a+ e+ V{(w@t+e)—4an}, 2 2. 2. (44) 
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the two infinite series arising from the alternative sign. For extreme 
values of 7 the roots are readily obtained by using the appropriate 
binomial expansion of the term under the square root sign in (4.4). 
Proceeding in this way with the aid of fig. 1 it is not difficult to show that 
ify is very large, that is the reaction is very rapid compared with diffusion, 
the terms in the general solution for M,/M,, which arise by taking the 
negative sign in (4.4) vanish, as does also the term from the imaginary 
root. The terms from the positive sign lead to 


MS 2a(1-++«) exp (—Bk,,?a*t) 
where the k,,a’s approach the roots of 
tan k,a=—«k,a , ae ates teeta dl (40) 
and 
B=D/(R+1)a?. eee meer cnt) 


This is the solution obtained by Wilson (1948) for the case of an 
immobilizing reaction which is rapid compared with diffusion. 

If, however, 7 is small because the reaction is infinitely slow (u=0) we 
find that the terms arising from the positive sign in (4.4) vanish provided 
a/l is not zero, which case is treated separately in§ 5. The terms from the 
negative sign combine with that from the imaginary root to give 


Le (1 Tot )-2 Se exp ae 5 ics 


M, \ R+1i+a i 1+Va-+ (I/a)*k,?a2 
where the k,,’s are given by 
tan k,a=—(l/a) k,a . Pe ne, (4,9) 


The whole term in the second bracket in (4.8) is to be recognized as the 
expression for simple diffusion from a finite bath i.e. diffusion in the 
absence of any immobilized component. Equation (4.8) also describes a 
simple diffusion process, therefore, and M,/M,, changes from zero at t=@ 
_ to 1—Rl/(R+1)(/+-a) at t= 00, which is easily shown to be the fractional 

uptake of solute to be expected in the absence of immobilized solute. 
Thus (4.8) indicates the behaviour to be expected on general argument, 
namely that for an infinitely slow reaction the sheet takes up, by simple 
diffusion, only the fraction of solute which it can accommodate in the 
freely diffusing state and none in the immobilized state. 

If on the other hand pa?/D is small because D is very large, all terms in 
the general solution vanish except the one associated with the imaginary 


root and we are left with 


M, fie Ra ) \ 

a == aa Le oe ot se) 4A) 

Ae. Raa? { (ee e cay) 
This expression is readily deduced from elementary considerations when 


diffusion is so rapid that the concentration of solute is effectively uniform 
through the sheet at all times. 
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The type of behaviour observed in a practical system for which pa*/D is 
very small depends on the time scale of the experiment. If this is such 
that the reaction occurs very slowly compared with the duration of the 
experiment the simple diffusion behaviour of (4.8) is observed. If on the 
other hand diffusion is very rapid compared with the time scale of the 
experiment the simple first-order reaction of (4.10) is observed. 


§5,. EvALUATION OF NUMERICAL VALUES 


When eqns. (3.9) and (3.10) are written in terms of p,/u and k,a we 
see that I/,/M., can be calculated as a function of Dt/a? if three parameters 
are known. The parameters are //a, that is the ratio of the volumes 
occupied by solution and sheet respectively, the partition factor R, and 
the modified rate constant for the reaction, wa?/D. Alternatively, since 
R=Ajp, a solution is defined by //a and the two rate constants pa?/D and 
da?/D. In some cases it is more useful to relate Rk to the fraction of the 
total amount of solute which is in the sheet finally, i.e. to M,,/ICp, by the © 
relations (3.11). For the cylinder the corresponding parameters are 
7a?/A, R and a?/D, and for the sphere 47a3/3V, R and pa?/D. 

Once the roots p,, k,a are obtained, the evaluation of each of the 
expressions for M,/M,, for the plane sheet, cylinder, and sphere is straight- 
forward provided /, A and V are finite. The cases of /, A and V infinite, 
however, need further consideration because the convergence of terms for 
which p, approaches —y can be very slow and numerical evaluation 
becomes awkward and laborious, particularly for small «a?/D. For the 
plane sheet when «= 00, (3.10) reduces to 


M, 2D*k,,2(1+ pyle)? EXP (Ppt) 
M. (R+1)p,207{(1+p,/u)2+R}? ~~. 


where now k,a=(n-+-3)7. As we saw in § 4 there are two infinite series in 
the general expression for M,/M.,. We shall confine attention for the 
moment to the series associated with the positive square root in (4.4) since 
these are the terms for which p, approaches — when Dk,,?a*/ya? is large. 
Substituting for 1--p,,/u from the second of (3.9) we find 


=1—2 (5.1) 
1 


My 5 PRAM Palle) OXP (Pal) 52) 
Me (RIO +p PP Rye (OM 
If p,=—p to the order of accuracy required, after the first » non-zero 
roots, we have ) 
Hs oy 5 RT Dall)® OXD LD aap eee ee ee 
M4 (RFI +p, PF Rika — RPI — 2px} exp (—ut 


; - (5.3) 
approximately since 


k,,2a* = 1 : e ° ° . . ° . (5.4) 
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The relationship (5.4) follows, for example, from (4.5) when «=o, 


since M,/M,,=0 whent=0. The error involved in use of the approximate 
form (5.3) is less than 


[See eXp (D440) Rk 5 
TEETER) TEP) |[ Fez]. . 6H 
and can be made as small as desired by choice of r. We may note in 
passing that since in (5.2) p,/u and 1+-p,/u occur only with R, this is a 
more convenient form of expression for computation than (5.1) when 
P,/e 18 small or near —1, particularly for large R, since it is less sensitive 
to the accuracy of the roots p,,/u. 

When l/a= «, it follows from (4.2 a) that there is no root for which 
k,a is imaginary. The complete expression for M,/M,, is therefore that of 
(5.3) together with terms arising from the negative sign in the roots of 
(4.4). On using (5.3) numerical evaluation of M,/M,, for l/a infinite is 
straightforward. The corresponding formulae for the plane sheet and the 
sphere are easily derived. 


§6. [IRREVERSIBLE REACTION 


A special case of the above solutions of particular interest is that of an 
irreversible reaction, when the rate of formation of immobilized solute is 
directly proportional to the concentration of free solute. In this case 
p»=0, but A is non-zero so that R= oo, «=0, M,,=IC,. The solution for 
the plane sheet, for example, for these values follows immediately from 
(3.9) and (3.10) provided J is finite. 

The solution for the case of /= oo is less obvious. When p»=0, the 
imaginary root (p,, k;) is given by 


lp;=— Dk; tanh k,a, k,?=(p,;+-AjD, . . ~~ (6.1) 
where k,=7k,’, and so when /= oo 
een Dea eee te te. act. (6.2) 


When «=0, »=0 and p, is small, we can expand exp (p;t) in powers of 
pt and write (3.10) as 


1IOo(1-+ pit IC'5 CXP (Pn 
Oe a 
: a Pi pla l ih a tele Pala’ 
l+o1—3pe8  2p% ‘To eDEs SDE at 9%, 
from which, when /= 00, so that lp, is given by (6.1) and (6.2), we have 
finally 


M, ODt 2.D?k,,” exp (Dat) 
Tone ard tanh g-+-4 sech? tb ee tanh 2 gees 


where 


a (6/4) 


k,a=(n-+4)0, b,2=—(Pytr/D, g=V(ra"/D). . « (6.5) 
The first term on the right-hand side of (6.4) gives the rate of uptake of 
solute due to the chemical reaction in the final steady state. 


820 J. Crank on the Simultaneous 


The forms of (3.12) and (3.14) for the cylinder, and of (3.16) and (3.17) 
for the sphere when p=0, «=0, are obvious. When A= © we have for 


the cylinder 


2 4D2k,2 Jt ‘ 
M, et ee yee exp (P ) (6.6) 
7*C a? *1(q) — 1*(9) Pn 
where . : 
Jo(k,jt)=0, &,2=—(Dar AD; q=v (Aa?/D).. a. Oa 
When V= ~ we have for the sphere 
3M 3Dt 3 6D?k,,? exp (ppt) 
ina, BoA (q coth gq—1)—$ cosech? g + Bj cule ar oshamenreres any ; 
(6.8) 
where 
k,a=nr, k,2=—(p_taA)/D, q=V (A/D). . « « (6.9) 


The solutions (6.4), (6.6), (6.8) can of course by obtained directly by use 
of the Laplace transform or otherwise. Danckwerts (1951) has shown 
that solutions for diffusion and simultaneous irreversible reaction can 
sometimes be obtained by a simple transformation of the expressions for 
the corresponding diffusion problem without reaction. 


§ 7, DESORPTION 
We have considered diffusion into a plane sheet initially free of solute. 
There is the complementary problem in which all the solute is initially 
uniformly distributed through the sheet and subsequently diffuses out 
into the solution. If free and immobilized solute are considered to be 
initially in equilibrium everywhere in the sheet, it is easily seen that the 
mathematical solutions presented above for absorption also describe 
desorption, provided M, is taken to mean the amount of solute leaving 
the sheet up to time t, and M,, the corresponding amount after infinite 
time. For desorption from a plane sheet, for example, we want solutions 
of the eqns. (2.1) and (2.2) satisfying the condition (2.4) but with (2.3) 
replaced by 
c=), 8=Sy, —a<e<a, t=0.. ye ees 
Writing 
Cy=Cy—c, 8,=S)—8, Pee ere ie) fF) 
it is easy to see that (7.1) and the other equations for desorption are 
identical with (2.1), (2.2), (2.3) and (2.4) with c,, s, written for c, s 
respectively, remembering only that AC)—S,)—0, if we have equilibrium 
throughout the sheet at t=0. Hence the equations and solutions for 
desorption are identical with those for absorption. 


§8. CALCULATED RESULTS 
Evaluation of the expressions of § 6 for the irreversible reaction and any 
particular set of parameters is comparatively simple and straightforward. 
On the other hand a considerable amount of painstaking labour is involved 
in the evaluation of the formulae for the reversible reaction even when 
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there is..an infinite amount of solute, and it is therefore worth while to 
present some numerical values. In tables 1, 2, 3 and 4 calculated values of 
M,/M,, are tabulated for the plane sheet, cylinder, and sphere. The 
arrangement of the tables was decided by ease of presentation. Three 
values of the partition factor R are included and for each R three rates of 
reaction. All tabulated solutions refer to an infinite amount of diffusing 
substance. The values of M,/M., are believed to be correct to within 
+1 in the third decimal place. 


, Table 1. Values of M,/M,, for na?/D—0-01 


plane sheet cylinder 


0-024 | 0-009] 0- : 0-057 
0-034] 0-013] 0-0: : 0-073 
0-048} 0-021] 0- 

0-067 | 0-036 

0-079 | 0-051 

0-086 | 0-065 

0-091 | 0-079 

0-099 | 0-114 

0-104} 0-147 

0-114} 0-210 

0-123] 0-268 

0-133 | 0-322 

0-178] 0-537 

0-220 

0-261 

0-299 

0-463 

0-588} 0- 

0-685] O- : 1-000 
0-815} 1- 

0-891 

0-987 

0-998 

1-000 


Solutions for the plane sheet are illustrated in figs. 2,° and4. Those for 
the cylinder and sphere show the same general features, differing only in 
detail. Figure 2 shows how the absorption curves change in shape and 
position as a?/D is varied between the two extremes given by = 00 
(infinitely rapid reaction) and »=0 (simple diffusion with no reaction). 
By plotting against (ut)! as in fig. 3 we can show the approach to the 
curve for D= © (purely reaction controlled) as j.a?/D tends to zero because 
of D becoming large. We see also in fig. 3 that the discontinuity in the 
gradient of the curve for D infinite appears as a ° shoulder * when 


Table 2. Values of M,/M,, for pa?/D=0-1 


plane sheet cylinder sphere 


R 1 10 100 1 10 100 1 10 100 


Dt/(R--1)a? 

0-005 0-057} 0-025} 0-017] 0-108] 0-044] 0-029} 0-155] 0-059 0-037 
0-01 0-080 | 0-035} 0-032] 0-149] 0-060} 0-054] 0-210] 0-076 | 0-068 
0-02 0-113] 0-052} 0-062] 0-205] 0-081] 0-101} 0-279] 0-095 | 0-127 
0-04 0-160] 0-076] 0-116] 0-277] 0-106} 0-188} 0-361] 0-117] 0-232 
0-06 0-197} 0-095] 0-166} 0-327] 0-125] 0-265 | 0-409] 0-136] 0-324 
0-08 0-227] 0-112] 0-211] 0-365] 0-143] 0-333 | 0-441] 0-153} 0-404 
0-10 0-254] 0-128] 0-253] 0-395] 0-159] 0-394] 0-463] 0-171] 0-474 
0-15 0-310} 0-164] 0:343| 0-445] 0-200} 0-519] 0-494] 0-212] 0-613 
0-2 0-354.| 0-198} 0-418] 0-476] 0-238] 0-615] 0-508] 0-252] 0-712 
0-3 0-418] 0-261} 0-536] 0-507] 0-309] 0-748] 0-522] 0-325] 0-838 | 
0-4 0-459} 0-320] 0-626] 0-524] 0-373] 0-832] 0-532] 0-391 | 0-907 
0-5 0-488} 0-374] 0-696] 0-535] 0-432] 0-887] 0-542] 0-451] 0-946 
1-0 0:559| 0-585] 0-889! 0-579] 0-652} 0-983] 0-585] 0-672] 0-996 
1-5 0-601 | 0-724] 0-959} 0-619] 0-786] 0-997 | 0-624] 0-804} 1-000 
2:0 0-638 | 0-816] 0-985] 0-655] 0-868} 1-000 | 0-660 | 0°883 | 
2-5 0-671 | 0-877] 0-994] 0-687] 0-919 0-692 | 0-930 
5:0 0-797] 0-983] 1-000} 0-809} 0-993 0-813] 0-995 
75 0-875 | 0-998 0-884} 0-999 0-886] 1-000 

10 0-923] 1-000 0-929} 1-000 0-931 

15 0-971 0-974 0-974 

20 0-989 0-990 0-990 


40 1-000 1-000 1-000 


Table 3. Values of M,/M,, for na?/D=1-0 


plane sheet cylinder 


R 1 10 100 1 10 100 1 10 100 


Dt/(R-+-1)a? 


0-005 0-057 | 0-028] 0-045 | 0-109] 0-051] 0-085] 0-155] 0-070] 0-120 
0-01 C-O81 | 0-045} 0-O8L] 0-151] 0-079] 0-151] 0-211] 0-105] 0-212 
0-02 0-115} 0-075] 0-136] 0-208] 0-129] 0-250} 0-283] 0-166] 0-346 
0-04 0-164] 0-131] 0-210] 0-285] 0-217] 0-378] 0-371] 0-273] 0-510 
0-06 0-203 | 0-181] 0-264] 0-341] 0-295] 0-465] 0-429! 0-366] 0-613 
0-08 0-237 | 0-226] 0-309] 0-385] 0-363] 0-531] 0-471] 0-445] 0-687 | 
0-10 0-267] 0-268 | 0-348 | 0-423] 0-423] 0-587] 0-504] 0-514] 0-743 | 
0-15 0-334 | 0-358] 0-480] 0-495] 0-547] 0-692] 0-564] 0-648] 0-840 
0-2 0-390] 0-432] 0-497] 0-549] 0-638] 0-767] 0-608] 0-743] 0-899 
0-3 0-481] 0-549] 0-607] 0-630] 0-768] 0-866] 0-676] 0-860] 0-959 
O-4 0-553] 0-637] 0-691] 0-691] 0-848] 0-923] 0-732} 0-922] 0-983 
O-5 0-613] 0-707] 0:758] 0-741] 0-898] 0-955 | 0-777 0-956} 0-993 
1-0 0-804] 0-895] 0-927] 0-891] 0-986] 0-997] 0-912] 0-997] 1-000 
1-5 0-900 | 0-962] 0-978] 0-954} 0-998] 1-000] 0-965] 1-000 

2-0 0-948 | 0-986} 0-994} 0-980} 1-000 0-986 

2-5 0-973] 0-993] 0-998} 0-992 0-995 

5-0 0-999] 1-000] 1-000] 1-000 1-000 

75 1-000 


Ene ens Ree rer el Sl 


0°75 


™,/ mM 


0:50 


0:25 


Table 4. Values of M,/M,, for pa?/D=10 


plane sheet | cylinder | sphere 


i 1 


Dt/(R+1)a2 


0-005 0-160 
0-01 0-224 
0-02 0-316 
- 0-448 
0-546 
0-623 
0-686 
0-796 
0-865 
0-939 
0-972 
0-987 
1-000 


SCL bo 
S Ol OS Or 


0-5 2:0 


7 F |: 
Coe(R+i)a3}! 


| ord. | My/ Mc. | abs. | {Dt/(R+-1)a2}? 


Plane sheet, R=10. Curves are labelled with values of pa?/D. 
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pa?/D=0-01. As a?/D is increased further the shoulder disappears ae 
an absorption curve with a point of inflexion. At still higher wa |D the 
inflexion becomes less noticeable as in fig. 2 and the curves have the simple 
shape commonly associated with diffusion. Figure 4 shows the rea 
of the parameter R, the partition factor between immobilized and ree 
solute. As R is increased the height of the shoulder decreases and if 
curves for R—100 were plotted on the present scale no shoulder would be 
detected for any value of a?/D. 


Fig. 3 


1-00 


0-75 


0-25 


0 0:5 lO Cutty "5 2:0 2:5 


MMe (ut) 


Plane sheet, R=10. Curves are labelled with values of pa?/D. 


One interesting feature of these results is that they indicate limits to the 
relative rates of diffusion and reaction outside of which the reaction is 
effectively infinitely rapid or infinitely slow as the case might be. Thus 
for the plane sheet the values of M,/M,, for .a?/D=10 differ by only a few 
per cent from those for an infinitely rapid reaction. The differences are 
greatest for R=1 where they are about 5%. At the other extreme the 
solution for na?/D—0-01, R=1 is the same as that for an infinitely slow 
reaction to the same degree of accuracy except at small times (see fig. 4 


for example). The differences increase in this case as R increases, being 
about 20°, for R=100. 
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The significance of these limits is perhaps easier to appreciate when they 
are expressed in terms of the half-times of the simple diffusion and simple 
reaction processes respectively. For simple diffusion into the plane sheet 
from an infinite amount of solution it is easy to show that the half-time, fj, 
that is the time at which M,/M,,=, is given by 

Dtqla?=—(4/n®) In w?2/16=0-2 . . . . . (8.1) 
approximately. On the other hand if immobilized solute is formed from a 
uniform constant, concentration, Cy, of free solute according to the 
equation 


Ga CieN ie mate ae fn 18.8) 
Fig. 4 
1-00 : 
0-75 
0:50 
M2 /Mo 
0-25 
O 0-5 ora. Te 2-0 2-5 


ord. | M,/Mu | abs. tye” 


Plane sheet, R=1. Curves are labelled with values of pa?/D. 


this reaction has a half-time, ¢,, given bv 

fies 2 0-Te eee ws (8.8) 
approximately, and combining (8.1) and (8.2) we have for the ratio of the 
half-tim es, 
fi 2S OD dy eam lay os. Comet (84) 
Thus we can say roughly that if the diffusion is more than a thousand times 
faster than the reaction, the behaviour of the diffusion-reaction process 
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is roughly the same as it would be if diffusion were infinitely rapid. On 
the other hand if the half-times for diffusion and reaction are comparable, 
the behaviour approximates to that for any infinitely rapid reaction. 
These statements indicate orders of magnitude only. They apply also to 
the cylinder and sphere. 

In conclusion the author wishes to acknowledge the computational 
assistance of Mrs. D. M. Lambert to whom the tabulated numerical values 
are due. 
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ABSTRACT 


It is shown that slight anomalies in the layer pattern near a dislocation 
can often be interpreted only as due to movement of the dislocation over 
a short distance. Examples are also shown in which a single dislocation 
has moved and cut across neighbouring layers. : 

Characteristic types of dislocation groups, termed limited slip-zones, 
which are often observed, are shown to be explicable, in all details, in 
terms of the ‘ generation ’ of dislocations by Frank—Read sources. The 
type of layer structure to be observed while such a source is still active, 
or has just ceased to be active, is deduced and some of the cases observed 
are shown. 


$1. INTRODUCTION 


Tue patterns of unimolecular layers observed on beryl crystals are static 
and necessarily display the situation extant when the layers finally come 
to rest, whether the final stage was growth or dissolution. The purpose 
of this paper is to present evidence that some of the dislocations present 
have moved in some stage of the history of the crystal ; this is deduced 
from the static picture through the observations of anomalies in the layer 
structure near the dislocations which are inexplicable in terms of either 
growth or dissolution ; alternatively, certain precisely defined types of 
dislocation groups are observed which give evidence for the existence 
of Frank—Read sources. Evidence is presented below to show that on 
the various crystals observed deformation has occurred, including both 
the movement of single dislocations over very small distances, and 
large-scale production of dislocations from sources in the interior, 
resulting in the displacement of considerable sections of the crystal. 
_. If any movement of a dislocation over a short distance is to be detected 
it-is ‘necessary to know 

(a) the dependence of growth rate upon orientation for the particular 
crystal face ; 

(b) the normal type of cross-linking which takes place between a 
dislocation and advancing layers. 

The examples cited in $2 were all observed on the prism faces of 
one beryl crystal. These faces had a rate of growth parallel to the c-axis 
of the crystal, which was usually some two to four times greater than the 
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rate of growth perpendicular to this axis. The normal type of. cross- 
linking with a single dislocation (Griffin 1951) is shown for a layer edge 
running perpendicular to the c-axis in figs. 1a and 1b, where the shading 
denotes the low side of the steps. 

The type of cross-linking is precisely similar for right- and left-handed 


screw dislocations. Propagation of layers past single dislocations is, in 
all cases, a simple repeat of stages (i) to (iv) in fig. 1. Inspection of this 


Fig. 1 


b 
(i) (ii) (iii) (iv) 
The normal type of cross-linking which takes place with a single dislocation. 
Both a right-handed and a left-handed screw dislocation are shown, the 
layer topography in the one case merely being a mirror image of the other. 


figure also shows that for the layer which passes on after cross-linking, 
the “ bump’ or kink in the edge accurately locates the previous position 
of the dislocation. 
A case of cross-linking which may appear identical to fig. 1 in its 
early stages is shown in figs. 2a and 2b. Here cross-linking results in 
(a) the formation of a depression which will be filled in rapidly ; 


(6) attachment of the dislocation to the type of layer edge shown 
in fig. 1. 
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Tt is obvious that layers of the stage (i) type of fig. 2 cannot persist 
after the first cross-linking, and that the second, and all subsequent 
cross-linkings, will be of the type shown in fig. 1, the dislocation pers 
attached to a normal type of layer. 

This can be generalized, and any complicated type of layer running 
from a dislocation can be considered. This is shown in fig. 3. 

We are led to the following general rule. No matter what the shape 
of the original layer attached to the dislocation, this layer will conform 
to the layer type of the incident system after the first cross-linking. This 
will also be true for all subsequent cross-linkings. 


Fig. 2 


(i) (ii) (iti) (iv) 


6 
(i) (ii) (iii) (iv) 

Stage (i) of both a and 6 shows a situation which under the microscope appears 
to be identical with that shown in fig. 1. However, as the steps attached 
to the dislocations are facing in the opposite direction in this case, the 
result of cross-linking is quite different. 


Alternatively this rule can be stated as follows: provided that the 
dislocation remains stationary, and that the conditions are not 
fluctuating markedly, all cross-linkings after the first will be of the types 
shown in fig. 1. 

We shall now consider some experimental observations, 


212 
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§ 2. MoveMENT OF A DISLOCATION OVER A SHORT DISTANCE 


Example 1 
This is shown in fig. 4, Plate LIV. Two such events occur on the 


plate ; they are isolated and shown in figs. 5a and 5c. (in all figures 
the c-axis is approximately horizontal unless otherwise specified.) 
In fig. 5a note (i) the re-entrant in the layer edge numbered 1, (ii) the 


position of the dislocation. 
Fig. 3 


An arbitrarily shaped layer is attached to the dislocation in the first instance, 
but the result of cross-linking is to attach the dislocation to a normal type 


of layer. 
Fig. 5- 
! le 
a b , 


a. This shows the layer shapes immediately adjacent to the lower dislocation 
which is indicated on fig. 4, Plate LIV, and will help to identify this area. 

4. The assumed layer shapes at the moment of cross-linking. The cross 
marks the position of the dislocation which is actually observed. 

ec. This shows the layer shapes close to the higher dislocation indicated in 
fig. 4, Plate LIV. Layer 1 is also common to fig. 5a. 


It is known that sharp re-entrants are eliminated at a speed considerably 
greater than the normal growth rate, probably for a re-entrant of this 
curvature about twice as fast. However, in the case of fig. 5 the direction 
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of elimination of the re-entrant corresponds to a minimum growth rate, 
and in this region the rate of growth parallel to the c-axis is some four times 
greater than the minimum growth rate. Consequently a re-entrant 
reaching the dislocation is unlikely to have been eliminated to the extent. 
apparent in this figure. 

The observed shapes of the layer edges can be satisfactorily explained 
if it is assumed that at the moment of cross-linking of edges 1 and 2 
the position was as shown in fig. 5b. Here the cross marks the observed 
position of the dislocation. It is shown above that such an edge-shape 
could not persist. Therefore this deviation from a normal shape could 
only have been produced between the penultimate and final cross-linkings. - 
of this dislocation before growth ceased. Thus the motion must have 
taken place just before growth ceased. 

In fig. 4, immediately above the layers just considered, is the situation. 
shown in fig. 5c. Applying similar arguments, it is obvious that layer O 
must have taken this outline just before growth ceased. This layer can 
have undergone no cross-linking. Note also that the bump on the lower 
part of layer O is entirely consistent with normal cross-linking of the 
dislocation, now attaching layer 2 of fig. 5a with layer O: this was 
the second cross-linking before growth ceased. Moreover, the position 
of this bump is such that the dislocation must have been near the top 
of the right-hand straight portion of layer 1 as seen in fig. 5b. The 
short straight sections of the layer edges running up to the dislocations 
in figs. 5b and 5c are not consistent with the normal growth rates. 
These normal growth rates are observed to hold for the layers immediately 
adjacent to those under consideration. 

The only explanation which can satisfactorily account for this 
behaviour is that the two dislocations have moved a short distance 
along a line perpendicular to the c-axis. Thus it appears that these two 
dislocations of the same hand which are 10 microns apart have each 
moved a distance of approximately 5 microns. 


Hxample 2 


This case is shown in fig. 6, Plate LIV. This could correspond to 
either fig. 7a or 76. In fact the shape of layer 1 shows that cross- 
linking took place just before growth ceased and, therefore, the right-hand. 
upper part of layer 1 and the right-hand part of layer 2 must correspond. 
Therefore the situation shown in fig. 76 must be that observed. 

The shape of layer 2 could not possibly be due to cross-linking, and, 
in fact, this shape would disappear at the first cross-linking. This shape 
is also entirely at variance with the observed growth rates and growth 
behaviour. The sharpness of the right-hand corner of the bump on 
layer 1 indicates that cross-linking can only just have taken place before 
growth ceased. Thus the only feasible explanation is a movement of the 
dislocation perpendicular to the c-axis, as shown in fig. 7c. The sharp 
corner resulting from this movement would soon be rounded off by 
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dissolution to the kind of shape observed. In this case the distance 
moved was probably some 6 microns, part of this movement possibly 
taking place after growth had ceased. 


Examples 3 and 4 

These examples are shown in figs. 8 and 9, Plate LIV. In this 
case interest is centred in the hook-shaped layer edges. In principle 
these examples are similar to example 2. In detail they differ in that the 
‘hook’ in these two examples is facing the advancing layer system. 
Both hook-shapes must have been produced since the last cross-linking 
took place. In fact in both cases the position of the bump on the last 
layer to be propagated past the dislocation suggests that the dislocation 
was located at the lowest point of the hook at the moment of cross-linking. 

Once again these examples can only be satisfactorily. explained if, 
since the last cross-linking, the dislocations have moved a distance 
perpendicular to the c-axis of 6 microns and 8 microns respectively. 


ON 


a and b The shading here denotes the low side of:the steps, and the layer shapes 
indicated in fig. 6, Plate LIV, could correspond to either of these two cases. 

c. The full lines represent the original and typical layer shapes, while the 
dotted line indicates the path of movement of the dislocation to its observed 
position at the point indicated by a cross. 


Some layers shown in fig. 9 are isolated in fig. 10. It is necessary to 
decide whether the dislocation is attached to layer 1 or layer 2. If the 
dislocation is attached to layer 1, it means that cross-linking has still to 
take place, but if the dislocation is attached to layer 2 then cross-linkin, 
must have taken place. : 

It appears that cross-linking has taken place because 

(a) the bump on edge 1 is the normal type produced by a little growth 
subsequent to cross-linking : bumps of this type do not result from 
normal growth of a layer attached to a dislocation ; 
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(b) layer 1 as far as the start of the bump, and layer 2 above this 
level appear to be co-linear. 

Further support for this view was obtained from observations with an 
objective of N.A.1.32. The dislocation and layer 1 were then almost 
resolved and appeared to be separate. 

The bump on layer 1 indicates the approximate position of the 
dislocation when cross-linking took place. At some time after cross- 
linking, therefore, the dislocation must have moved 10 microns to its 
observed position. The layer edge running from the assumed position 
of the dislocation at the moment of cross-linking to the position observed 
in fig. 9 is only perpendicular to the axis in its lower part. This could 
be due to a small amount of growth taking place after the dislocation 
had moved an initial short distance perpendicular to the c-axis, further 
movement in the same direction taking place after growth had ceased. 
It is probably coincidental that the dislocation has come to rest when 
almost in contact with layer 1. 


Fig. 10 


Z 


The observed layer shapes in the area marked by the left-hand pair of vertical 
and horizontal arrows in fig. 9, Plate LIV. 


Example 5 

The layers indicated in fig. 11, Plate LV, show three ‘ hooks’ 
which are very similar to the examples shown in figs. 8 and 9. The 
layer system is advancing from left to right, so that the * hooks’ are 
again facing the incident layers. Pu 

Again these cases can only be satisfactorily explained if the three 
dislocations which are all of the same hand, have moved distances of 
approximately 6 microns, all in the same direction. 


EHzample 6 

This is shown in fig. 12, Plate LV. The events leading to the 
observed layer shapes can be represented as in fig. 13. Once again the 
only satisfactory explanation of these shapes is the type of movement 
shown in the figure. Figure 13 (i) shows in full line what may be assumed 
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to be the original shape of the layers, while the dotted lines show the path 
of movement of the dislocations. Very slight growth subsequent to such 
movement will produce the layer shapes shown in fig. 13 (ii), these being 
identical with the observed shapes. This movement over a distance of 
2 microns and 3 microns respectively is normal to the c-axis, as usual, 
In this case movement probably occurred only a very short time before 
growth ceased. 


Example 7 
This is shown in fig. 14, Plate LV, and the processes involved are shown 
in fig. 15. (This latter figure may be more readily comprehensible if it is 
remembered that although the eye moves from left to right for the 
successive stages the direction of advance of the layer system is from 
right to left.) 
Fig. 13 


- 


€--- 


(i) (i) 
The area of interest in fig. 12, Plate LV, is shown on the right. On the left 
the full lines represent assumed typical layer shapes, while the dotted 
lines represent the path of movement of the dislocations. Slight growth 


on the new layer edges so produced would th F 
shown on the right. F ou en result in layer shapes as 


It is assumed that the original situation was as shown in fig. 15 (i), 
and that the two dislocations moved just in front of the incident yee 
causing the latter to retreat a little (e.g. because of the heat Sper 
in dislocation motion) as indicated in figs. 15 (ii) and (iii). Eventually 
after an abnormally large kink had been formed, cross-linking took place, 
producing the final configuration shown in fig. 15 (iv), which corresponds 
to that observed. The layer edges running from the main bumps down to. 
the dislocations must represent the paths of the latter. These paths make 


an angle with the main part of the layer edge and 
perpendicular to the c-axis, SE Ldaeth 23 ee 
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Hxample 8 


This is shown in fig. 16, Plate LV. The two twin-dislocation 
loops indicated are of quite anomalous shapes. These shapes are totally 
inconsistent with the observed growth rates on this face, and can only be 
satisfactorily explained in terms of movement of the lower dislocation of 
each loop. 

Example 9 


See fig. 17, Plate LVI. This area is very close to that shown 
in fig. 16. There are two layer edges, indicated in the centre of the figure, 
which have short straight sections before the dislocation. The shape of 
these sections is inconsistent with the generally applicable growth rates 
and can only be explained, once again, by movement of the dislocations 
perpendicular to the c-axis. This movement probably took place after 


Fig. 15 


I 


(i) (ii) (iii) (iv) 

Tn this figure although the successive stages are drawn from left to right the 
layers are advancing from right to left. Stage (iv) shows layer shapes 
corresponding to those indicated in fig. 14, Plate LV. The situation shown 
in stage (i) is entirely typical, but in (ii) and (iii) the moving dislocations 
cause an ever-increasing length of kink, until cross-linking produces the 
layer shapes shown in (iv). 


growth had ceased. The dislocation at the termination of the upper 
arrowed layer edge must also have moved upwards. The way in which 
further growth would have resulted in these layer edges acquiring normal 
shapes has already been fully discussed (Griffin 1951). 


Example 10 

In this example, shown in fig. 18, Plate LVI, all three of the 
dislocations shown appear to have moved. The previous positions of 
the dislocations can be located from the ‘ bumps’ on the last edges te 
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have undergone cross-linking. In each case the distance moved 1s about 
4 microns. The sections of the layer edges marking this movement 
are approximately perpendicular to the c-axis, the deviation from true 
perpendicularity probably being due to a slight amount of growth taking 
place after the dislocations had moved. 


General 

The examples given above comprise only some of the cases of this 
type of movement which have been observed. Movement along a line 
perpendicular to the c-axis of the crystal corresponds to the trace of a 
(0001) plane in the face, and this plane may well be one of the 
characteristic slip planes of the crystal. The other family of slip planes 
which one might expect are the first order prism faces {1010}. These 
planes would, of course, have a trace parallel to the c-axis in this case. 


Fig. 20 


2 


The area of fig. 19, Plate LVII, which is of interest. 


The straight edge produced by such a movement is parallel to a 
close-packed direction (i.e. it is perpendicular to a direction of minimum 
growth rate) and the types of edge produced by slight dislocation 
movements and by normal growth processes are normally almost 
indistinguishable. It is therefore difficult to find examples which can be 
ascribed with certainty to dislocation movement. However, two cases 
have been observed in which, because of the local topography, it appeared 
probable that movement parallel to the c-axis had taken place. 


§ 3. MovEMENT oF a DISLOCATION ACROSS ONE OR MORE LAYER EDGES 
Example 11 


See fig. 19, Plate LVI. In this case interest is centred on the two 
layer edges isolated in fig. 20. Note that: (a) the right-hand part of 
layer 2 continues across a patch which partially obscures it ; there is no 
doubt, however, that layer 2 is attached to the dislocation as shown in the 
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figure ; (6) that section AB of layer 2, when produced, very nearly touches 
the lowest point of edge 1. The observed layer shapes can then be 
explained in terms of the processes shown in (i) to (iii) of fig. 21. 

In fig. 21 (i) layer edges 1 and 2 are of typical shape. Movement of 
the dislocation attached to layer 2, as shown in (ii), will produce the 
situation shown in (iii). It should be noted that because the steps in (ii) 
appear to cross this does not mean that the height of the steps is 
multi-valued at this point. The steps are, in fact, still distinct, although 
not resolved at Q. This can be explained as follows. 

At stage (i) of fig. 21 the levels on opposite sides of the steps differ 
in height by one unit cell, and the layer enclosed by edge | is one unit 
cell higher than the higher ‘level bounded by layer edge 2. If the 


Fig. 21 
- 
| 
S) 
R 
Q 


‘a : e 
(i) | (ii) (iii) 
Starting from typical spiral layer shapes, as shown in (i), movement of the lower 


dislocation, as in (ii), can produce the layer shapes shown in (ili), which 
correspond closely to those observed in fig. 19, Plate LVIT. 


dislocation of layer 2 moves towards layer edge 1, as in (ii), then the 
increasing length of edge attached to the dislocation still marks a difference 
in height of one molecule. This applies to an area in which, previous to 
the movement of the dislocation, no such difference in height had existed. 
The movement of the dislocation, in fact, has caused a depression of the 
lattice on one side of the step. As soon as the dislocation crosses edge 1 
the lattice to the right of the dislocation, and hence the step QR, is 
depressed by one lattice spacing. Therefore PQ of layer 2 and QR of 
layer 1 (see (ii)) are at the same level, mark steps of the same height, and 
meet at a point. Similarly, TQ of layer 2 and QS of layer | are now at 
the same level, mark steps of the same height, and meet at a point. In 
other words, PQR and TQS are now distinct steps which touch at a 


point Q. 
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If, after the situation shown in (ii) has been reached, a slight amount of 
growth or dissolution takes place, then the two layer edges will move 
apart. Dissolution will take place at these very sharp corners at Q, so 
as to round these corners off to a curvature not greater than that of the 
‘critical nucleus, and this, followed by a slight amount of growth, would 
produce the situation shown in (iii), which can be compared with fig. 19. 

It is possible for layer shapes of the type of 1 and 2 to be produced 
by dissolution alone. Such a mechanism cannot have been responsible 
in this case because of (a) slight, but important, differences of edge shape 
very close to the dislocations, (b) the closeness of other layers, which show 
no signs of dissolution themselves, and which would also prevent the 
ideal conditions necessary for the production by dissolution of layer 
shapes of this type. 


Fig. 22 


(i) (ii) (iii) 
The type of layer shapes produced by movement of a single dislocation across a 
number of layers. 


The distance moved by the dislocation was approximately 24 microns, 
The direction of movement has been partially masked by sight subsequent 
growth, but was almost perpendicular to the c-axis. 


Hxample 12 


Consider first the movement of a dislocation across a number of layers 
(fig. 22). This is merely an extension of the case considered immediately 
above. If the original position be as in (i) of fig. 22, with the levels of 
the layers as marked, movement of the dislocation downwards will give 
(ii), and subsequent slight dissolution at the sharp corners will 
produce (iii). Thus such circumstances will typically produce 
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(a) a series of kinked layer edges in which the short sections of the 
kinked edge are roughly co-linear for all the layers, and 

(6) layer-edge shapes such that the right- and left-hand sides of any 
kinked edge are respectively co-linear with one of the branches of either 
the upper or lower neighbouring kinked edges. 

These two conditions will become less obviously applicable as growth 
proceeds after movement of the dislocation, as the abrupt kinking will 
tend to be grown out, while the co-linearity of the now different layers 
will also gradually be destroyed. 

A situation almost identical with the single case of fig. 22 is shown in 
fig. 23, Plate LVII. Unfortunately the visibility is low in the area 
of interest which is contained by the arrows. Within this area the 
situation corresponds very closely to fig. 22 (iii), except that the layers 
have a decreasing spacing. One point of note is that the short kinks, 
which are strictly co-linear, lie perpendicular to the c-axis, and would 
normally be rapidly eliminated by growth. The dislocation therefore 
probably moved after growth had ceased. Although the significant 
features are of too poor visibility to show well in reproduction, it is 
possible to decide from the layer shapes that the dislocation probably 
moved upwards in fig. 16. The point where the dislocation came to rest 
cannot be decided with any certainty. 


§ 4. Limrrep Sure-ZonEs 


Some examples of limited slip-zones have already been shown 
(Griffin 1951). Another typical example is shown in fig. 24, Plate LVIT. 
As observed on a crystal surface these features consist of a line of 
dislocations which is divided about the centre into two equal groups of 
opposite hand. 

Limited slip-zones may conceivably be formed from the perfect lattice 
by slip taking place over part of a crystal plane. If the initial three- 
dimensional nucleus of a crystal is composed of perfect lattice it may be 
that processes of this kind result in the formation of dislocations in the 
lattice, thereby enabling the seed to continue to grow steadily at low 
supersaturations. However, the stresses necessary to produce slip in a 
perfect lattice are very high, and limited slip-zones are sufficiently common 
on some crystals of beryl for it to be obvious that most of them must 
have been formed when the lattice was no longer perfect, by some 
mechanism requiring much smaller stresses for the creation of dislocations. 
A mechanism admirably suited to account for all aspects of these features 
is provided by the multiplication process of Frank and Read (1950). 

In this process it is assumed that a dislocation line is anchored at 
two points. Under a sufficient stress this line will move outwards in the 
slip plane to encompass a larger area, and in doing so will eventually touch 
itself and so create a series of closed dislocation loops, in a way precisely 
similar to the behaviour of a double dislocation loop in the dislocation 
theory of growth. In this way a large number of dislocations can be 
generated from a single dislocation line. 


840 L. J. Griffin on the 


The shape and spacing of the loops so produced will be dependent on 
the stress, the elastic constants of the material, and the length of 
dislocation line which is acting as a source. If we consider a number of 
loops spreading out from a source and expanding in the slip-plane it is 
obvious that we have precisely the conditions for limited slip-zone 
formation. The dislocation loops will normally be compound in form 
(i.e. neither true edge nor true screw), so that intersection of a loop and 
a face will result in two points of emergence of dislocations with screw 
components of opposite hand (except in special cases). The two 
‘ dislocations ’ (or, strictly speaking, the two points of emergence of the 
same dislocation line) will continue to move outwards from the original 
point of contact while the ‘ source’ continues to be active. (Stages in 
this development are shown diagrammatically in fig. 25.) When the 


Fig. 25 


— — + + 


A section in the slip plane representing the action of a Frank—Read source in an 
isotropic material. The central arc represents the expanding ‘ source ’ 
dislocation line which produces the closed dislocation loops shown. Each 
of these loops expands to cut the crystal face in two points, the screw com- 
poner at the two points of intersection of each loop being of opposite 

and. 


dislocation lines come to a halt, either because of removal of the stress 
or other reasons, the ideal situation to be seen will consist of a movement 
step one unit high at the left- and right-hand extremities and steadily 
increasing in height towards the centre. 

The presence on the face of a closely spaced layer system before the 
source was active can result in a more complex situation. The various 
stages which would be observed on the face of the development of the 
slip-zone are shown diagrammatically in fig. 26. (Stages (i) to (iv) of 
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this figure are drawn, for the sake of clarity, as if dissolution was taking 
_ place simultaneously with the movement of the dislocations. This 
allows the connections between different layers to be shown, although 
in practice it would be impossible to resolve them.) Figure 26 (i) shows 
the situation when the first loop has just intersected the surface and the 
two dislocations are moving away from the point of contact. In fig. 26 (ii) 
the dislocations have cut the first two-layer edges, producing unit kinking. 
The original pair of dislocations are continuing to move outwards in 
fig. 26 (iii), leaving a series of layers behind them showing unit kinking, 
and meanwhile another dislocation loop has intersected the surface, the 


Fig. 26 


iter 


(i) (iii) 


ithe 


(iv) (v) 


The development of a limited slip-zone as observed on a crystal face with a layer 
system already present. 


two dislocations attaching edges as shown. These two dislocations are 
just about to cut other layer edges, and so produce two doubly kinked 
‘layers. In fig. 26 (iv) the slip-zone is limited on left and right by 
dislocations which are moving outward and contain @ series of layers 
showing unit kinking until the second pair of dislocations are met with. 
Between this pair all layers show double kinking as far as the centre of the 
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system, where the third pair of dislocations is now moving outwards and 

the left-hand dislocation of the pair is about to cut a layer edge and so 
produce triple kinking. (The asymmetry of attachment of this central 

dislocation pair is due to the asymmetrical nature of the kinking process.) 

The next pair of dislocations to be produced will be attached, in the same 

way, toa singly kinked and a triply kinked layer edge, producing quadruple 

kinking on intersection with other layer edges and so on. The type of 

layer structure which would normally be seen in practice after processes 

of this kind is shown in fig. 26 (v), where the plus and minus signs indicate 

the presence of dislocations. 

Further growth on this face will result in the advance of the main 
layer system, so that layers not attached to dislocations would have the 
kinks grown out by cross-linking and growth together, these layers 
eventually passing through the dislocation group, while at the same time 
linking of dislocations of opposite hand will take place (in the way 
discussed in a previous paper). Eventually, a layer system will be centred 
on the dislocation group as shown, for example, in fig. 24. (It is, of course, 
possible for the dislocation group to be completely “ dominated ° by the 
incident layer system, but this is not normally observed.) 

Thus it can be seen that the examples of limited slip-zones so far shown 
are derived from sources which were active before growth ceased. The 
processes shown in fig. 26 can only be observed if sources have been 
active either just before, or after, growth ceased. Examples of this 
kind are shown in figs. 27 and 28, Plate LVIII and fig. 29, Plate LVII. 

In fig. 27 the source appears to have been active along the line 
indicated, and in the original it is possible to trace the kinking of some of 
the layers. A layer count above and below the zone from the last 
dislocation on the left, to show definite evidence of layer kinking, gives a 
difference in count of two layers, although it is possible to distinguish 
a number of dislocations. 

The same type of behaviour is shown in fig. 28. On the right of 
the light patch which blocks out part of the zone it is again possible to 
trace the kinking between layers. (Unfortunately, in the general area 
in which figs. 27, 28 and 29 were observed the visibility of the layer edges 
was poor. Largely because of this no additional information could be 
obtained from observation with a 2mm oil immersion objective, as 
compared with the 4 mm objective of N.A. 0.95 with which these pictures 
were taken.) A layer count of this zone gives an equal number of layers 
above and below, so that the development is presumably of the ideal kind. 

A well-developed slip-zone is shown in fig. 29. At either end of the 
zone, where the degree of kinking is less, it is again possible to trace out 
the kinking between layers, and some dislocations can be seen. 

In all three cases the zones are approximately parallel to the c-axis of 
the crystal. Other features similar to these have been observed, but 
only in these cases was it possible to trace the connection between some 
of the layers in the zone. 
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On another crystal features were observed which may correspond to 
hmited slip-zones of considerable extent, in some cases extending from 
edge to edge of the face, and having been developed either just before, 
or after, the cessation of growth. The orientation of these features 
corresponded approximately to the trace of pyramidal planes in the face. 
However, there were always certain aspects of these features which 
could not be explained purely in terms of the movement of dislocations, 
so that similar evidence from other crystals is needed before any definite 
decision can be made about them. 

The number of dislocations present in limited slip-zones would be 
expected to vary widely, and this is found to be the case. The spacing 
of the dislocations will be largely governed by the factors influencing 
the shape and spacing of the dislocation groups, and is found to vary 
considerably. 

The distance apart of the two dislocation groups will also depend on the 
distance of the source from the face. If the source lies very close to the 
erystal face (and the other conditions are appropriate) the spacing 
between the two groups will be very small. On the other hand, if the 
source is far from the face the loop will be large by the time that it touches 
the face, so that a considerable section of the loop may be almost 
tangential to the face. This will result in the two groups being widely 
separated even if the original loops were closely spaced. When sufficiently 
far apart each group will be identical with, and act as, a lineage boundary. 
The table below gives approximate figures of the number of dislocations 
and their average spacing between dislocation groups for a number of the 
smaller limited slip-zones which have been observed. 


reir bar of Average spacing | Distance between Orientation 
: : of dislocations two groups relative to 
ee ous (microns) (microns) c-axis 
12 Led 12 parallel 
12 4-0) 15 parallel 
22 1-6 8 parallel 
24 4-] 20 parallel 
34 1-8 1-5 perpendicular 
64 1-5 27 perpendicular 
100 0-9 50 parallel 


Features which have long closely packed lines of dislocations, almost 
certainly in the form of limited slip-zones, are not uncommon, but these 
are not listed in the table as it is difficult to obtain accurate information 
from them. These features are marked by lines of * rods ’, a short stubby 
type of projection which will be discussed fully in a later paper in this 
series. For our present purpose it is only necessary to know that these 
rods only occur in areas which have a steep gradient parallel to the c-axis. 
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The size of the rods in any one area is directly related to the steepness 
of the gradient, or, in this case, the difference in height at a disconaae 
in level. These features are typically observed as a long line of “rods: 
producing a sudden break in the layer system, and in the simplest cases 
the rods increase steadily in size towards the centre of the feature, tailing 
off gradually towards each end. This indicates that the sudden change of 
level is greatest at the centre, and this can often be confirmed by 
examination of the layer topography. Many dislocations can often be 
observed lying along the bases of the ‘ rods ’, but it is often not possible 
to determine the sign of these dislocations. However, after the line of 
‘rods’ ceases the line of the feature is sometimes continued by a row 
of dislocations, all of the same hand. Thus the evidence is entirely in 
favour of these features consisting of either large single slip-zones or, 
in some cases, of series of contiguous slip-zones resulting from the 
simultaneous activity of several sources in the same, or neighbouring, 
planes. 

Features of this kind are too large for reproduction if the layers are 
to be satisfactorily resolved, but a feature which is in some ways similar 
is shown in fig. 30, Plate LIX. In this figure the magnification 
is too low for the layers to be satisfactorily resolved, but the main 
features can be distinguished. 

In this case two boundaries contain a strip, these boundaries being 
marked over most oftheir length by ‘rod’ structures which eventually 
end at the top of the figure in a more easily observed layer system. The 
right hand of the two boundaries starts first at the bottom of the figure 
and is distinguished by a line of small rods. It can be observed from the 
general nature of the layer structure that there is a steady decrease in 
height from the bottom to the top of the figure. Thus the sudden 
appearance of a line of rods can only be due to the commencement of a 
discontinuity in level due to the presence of a number of dislocations. 
(At higher magnifications it is possible to see that below this line of rods 
the boundary is marked by a line of thirty or more dislocations, all of the 
same hand.) When the right-hand boundary is roughly level with the 
start of the left-hand boundary the rods decrease in -size and then 
disappear, indicating an approximate equality in level between this 
boundary and the layer system to the right. The boundary is blank for a 
considerable length and the rods then again increase to a maximum size 
and then decrease again, indicating equivalent variations in level. 
After another short blank stretch the rods again appear and build up 
rapidly to become very large, showing that in this region the central 
strip must be much higher than the layer system on the right. 

At the lower end of the left-hand boundary the rods again point to 
the right, indicating that the layer system on the left is above the level 
of the strip. As one moves upwards this rod structure increases rapidly 
in size until it soon fills the whole of the central strip. The rods then 
decrease and eventually become very small before a blank section of the 
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boundary is reached. Another rod system then increases and decreases 
in size until another blank section of the boundary is reached. Above this 
point a rod structure commences which points to the left and gradually 
increases in size. This means that the central strip is now above the level 
of the layer system on the left (as it is also on the right). 

Figure 31 shows diagrammatically cross-sections of this area along the 
lines indicated in fig. 30. It can be seen that these variations in level 
of the strip cannot be correlated with layer spacing inside and outside 
the strip, as the layer spacing inside the strip is consistently larger and 
bears no relation to the layer systems on either side. Such variations in 
level can therefore only be due to the presence of groups of dislocations, 
and it seems probable that some form of slip process (presumably due to 


Approximate cross-sections across the feature shown in fig. 30, Plate Xe 


the simultaneous action of a number of different sources) has produced 
this curious ‘ buckled’ strip. Although the boundaries of the strip are 
mainly approximately perpendicular to the c-axis, there is one section 
where the deviation is considerable. 


§ 5, CONCLUSIONS 


The many slight movements of single dislocations which have been 
observed do not show, except over very small areas, any correlation 
between the sign of the dislocation and the direction of movement. As 
the energy of a dislocation is high, probably being several electron volts 
per atomic layer in beryl, there will be a considerable effective line 
tension on each dislocation. Consequently many of’ these slight 
movements probably correspond to a straightening and shortening of the 
dislocation lines, the movement possibly being aided by local internal 
stresses. 


3M2 
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The various stages, and degrees, of development of limited slip-zones 
which have been observed provide evidence in favour of the action of 
Frank-Read ‘ sources’ which is almost as direct and conclusive as one 
could expect to obtain from static surface studies. As dislocation 
‘sources’ of this type produce considerable amounts of slip in a material 
as strongly bound as beryl, it appears probable that this process is of 
application to a wide variety of substances. 

Particularly in the case of the movement of single dislocations, it is 
possible in some cases to fix a time scale, for, as has been shown, the 
movement must have taken place in the time interval between two 
cross-linkings. Unfortunately the rate of advance of a layer, which 
necessarily governs this time scale, must be assumed, and could vary 
between very wide limits, so that the time scale for any one crystal must 
remain entirely arbitrary. 

The basal plane (0001) and the first-order prism planes {1010} appear 
to be the common slip planes in beryl. Observations on (0001) faces of 
beryl have provided independent evidence of the {1010} family of slip 
planes. 
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ABSTRACT 


The possibility is considered that cosmic rays are accelerated in 
regions of the galaxy of intermediate particle density (+10 cm-%). 
It is found that a conservative estimate of the magnetic fields in such 
regions would lead to a radio noise flux of the observed order of magnitude, 
and that the observed spectrum could easily be produced. A possible 
cause is suggested for the low energy cut-off in the cosmic-ray spectrum. 
If the proposed mechanism accounts for any large proportion of the 
cosmic rays, these should be accompanied by a small flux of y quanta 
which should show the same anisotropy as the galactic radio noise. 


THE observed energy flux falling on the earth’s atmosphere as galactic 
radio noise of frequency between 2107 and 108 c/sec is of the order of 
3x 10-7 that of cosmic rays. This disparity has led several authors to 
seek to explain radio noise as an interaction of cosmic-ray particles with 
the galaxy. Radio frequency radiation due to radiative collisions 
between protons and electrons in interstellar gas gives much too weak 
a flux, and that from the Cherenkov effect of relativistic particles passing 
through hydrogen could only account for the observed radio flux under 
extremely improbable conditions of low temperature (<0-01°K) and 
ion density (proportion of ionized atom <10~*). Ginsburg (1951) has 
shown that the radio noise could only be produced by the interaction of 
cosmic rays with a general galactic magnetic field of 10-° gauss, if the 
flux of electrons of energy = 10° ev in the galaxy were some three times 
the proton flux observed in the primary cosmic rays. However, it has 
now been shown by Gold (1952) that the polarization of starlight, which 
was thought to be the strongest experimental evidence for the existence 
of such galactic fields, can be explained by other means. It is therefore 
uncertain whether such general fields exist in the galaxy, and it is of 
interest to examine the interaction between cosmic rays and regions of 
local field concentration. Ginsburg (1951) and Alfven and Herlofson 
(1950) have concluded that such regions would have to be assumed very 
different from our own solar system ; but this is to be expected, since the 
observed radio noise requires the existence of sources very much stronger 
than the sun. 

The most favourable assumption is that such regions are themselves 
the source of cosmic rays. In want of a non-electromagnetic theory of 
cosmic-ray acceleration, this seems indeed very probable. Kwal (1951) 
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has considered the case of a flux of high energy protons in a stellar envelope 
with a field of the order 10° gauss, and has concluded that this might 
explain the solar contribution to the radio noise. If a comparable flux 
of fast electrons were supposed to accompany the protons, the energy 
flux of galactic radio noise could also be explained in this way; but 
there would remain the difficulties : 

(a) that the observed radio sources do not coincide with bright visible 
objects ; and 

(b) that the ionization density in such regions would almost certainly 
be high enough to give rise to a strong circular polarization (cf. Ryle 1950) 
which is not observed in the galactic radio noise. 

It seems at least possible that cosmic-ray acceleration would take place 
in regions of comparatively low density such, for instance, as might be 
associated with nova shells or with condensations of matter in the early 
stages of star formation. It is of special interest to consider the 
production of radio noise in such regions since the identification by 
Bolton, Stanley and Slee (1949) of a radio source with the seat of a super 
nova explosion, the Crab nebula. We shall therefore consider what 
properties must be postulated for the regions so that the production in 
them of the observed cosmic-ray flux would simultaneously give rise to 
the galactic radio noise. , 


§ 1. PostuLaTED SouRcE REGIONS 


If the local magnetic fields were in thermodynamic equilibrium with the 
thermal motions of the matter, they might be large ; but in the absence 
of any complete theory of such fields we shall make the much more 
stringent assumption that they are so restricted that the radius of 
gyration of electrons in thermal motion is large compared to their mean 
free path for collisions involving the transfer of a large fraction of their 
energy (cf. Hoyle 1951). This leads to the approximate relation 


H <4nce’m1(3kT)-32N, 
where 7 is the temperature, NV the particle density and e and m the 
charge and mass of an electron. Thus if 7 is expressed in °K, then to 


order of magnitude 
we —12 fh 3/2 
H S Ome N ( 10! ) . . . . . . . ( 1 ) 


Most of the observed radio noise sources subtend angles less than 
6-8’ at the earth (though exceptions do occur (Bracewell 1952)), and at 
least some are at distances not less than 1018 cm (Smith 1951). If many 
of them are within our galaxy, they are presumably at distances 
comparable with interstellar distances, and may thus have linear 
dimensions D<101°—10!6em. For charged particles to be accelerated 
over a long path within the regions we require a mean free path for the 
collision of energetic particles, AX .D, and a radius of gyration R< D. 
These conditions could be satisfied, for particles of energy 10°<H<10" ey, 
by putting H>10-° ¢ and N<10" em-3, 
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«leaving unspecified the precise mechanism of particle acceleration, 
but only assuming it to be electromagnetic in nature, we must assume 
that the region would contain a flux of both protons and electrons of 
high energy. These would only be strictly bound within the region on 
rather special assumptions, but the tendency for charged particles to 
escape might well be small. The predominant mode oe escape would 
then be as the uncharged products of collisions—neutrons and y quanta. 
The quanta would subsequently nearly all escape from the galaxy without 
interaction, but those neutrons with H<101%ev would have a proper 
time for traversing the thickness of the galaxy long compared to their 
half-life for 6 decay, and would thus give rise to a cosmic-ray proton flux. 
Neutrons of energy H=10** ev would travel only about 10™ cm during 
a proper time equal to their half-life and would thus seldom escape from 
the accelerating regions. These assumptions would therefore lead 
naturally to the observed low energy cut-off in the cosmic-ray spectrum 
by a mechanism entirely different from that proposed by Unsold (1951), 
which would produce a flattening of the energy spectrum, rather than a 
cut-off, at low energies. The decay electrons would in general receive 
only a fraction of the neutron energy of the order of the ratio of the 
masses of electron and proton, and would therefore not contribute 
appreciably to the observed cosmic-ray flux. 

Since it is known (Perlow and Kissinger 1951) that quanta (at least 
those of low energy) form not more than a few per cent of the cosmic-ray 
flux, we must suppose that the mean path of a charged particle in the 
galaxy is much longer than that of a y quantum. We must, however, 
assume in any case that the proton paths are very much curved in the 
galaxy to account for the observed high degree of isotropy in cosmic rays. 
~ Since most of the cosmic-ray energy flux is carried by particles of 
H<10115 ev, and the radius of gyration of a particle of this energy in 
a field of 10-? gauss is R107!cm, which is comparable with the 
thickness of the galaxy, we need not assume any galactic field in excess 
of that required by relation (1) in order to ensure this for all but the 
higher energy particles. Ifthe most energetic cosmic-ray particles have a 
galactic origin, it is difficult to explain their isotropy without a large 
general galactic field. 


§2. Tue InTwRACTION OF ELECTRONS WITH THE FIELD 


The energy of the electron flux will be partly radiated into the galaxy 
as y quanta, of energy comparable to the electron energies, due to 
collisions with nuclei, and to a less extent as low energy quanta due to 
interaction with the field, some of which will fall in the radio noise 
frequencies. The latter process is described by the theory developed by 
Vladimirskii (1948) and Schwinger (1949). Thus the rate of energy 
radiation along the path of the electron is 


| | dE 2/2 \3/ B\ | = 
— ——— —. ———s 2 ° ° e ° ° . Pe 
(° dx (a) (<2) i’, (2) 
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where H is the component of the magnetic field perpendicular to the path 
(which we shall assume to be of the same order as |H|). Most of the 
radiation is at frequencies of order 
leH (ELE \? 
a (ee 
"1 Oo me (=a) (3) 
and the spectrum is roughly as shown in fig. 1 (which is plotted from figures 
given by Vladimirskii (1948)), the intensity varying with (v/v,)}* for v< 4. 
The ratio of the total low-frequency radiation flux to the y flux from 
electrons of energy # is 


WH). ab Nn 2 fren ee an 
Wo) 7 wes ne) 
Fig. 1 
P(y) 


%, 

The spectrum (p(v)) of electromagnetic radiation emitted by an electron 
moving in a magnetic field. 
where A is the mean free path for energy loss by radiative collision. 
The fraction of W(H), [W(R)]/[W()], which lies in the frequencies. 
2x 10’<v<108, can be obtained by integrating fig. 1, and from this we 
can find f(H)=[W(R)]/[W(y)]. For a fixed value of »,, [W(R)]/[W(A)] 
and #?H are constant, and, using eqn. (4), we can thus obtain a general 
relation between f and £ for all values of H and A by plotting fA-1H~?” 
against H'*H, ‘This is shown in fig. 2, in which the units are cm, gauss. 
and electron volts. 

The observed ratio of flux of radio noise to cosmic-ray y component 
will correspond to the value of f averaged over the energy spectrum of the 
electron flux within the sources. This spectrum is not known, for it 
cannot be assumed to be the same as the observed spectrum of primary 
cosmic rays, but it is plausible to assume that it is a steeply decreasing 
function of H. If, then, we consider only those electrons of such energy 
that Ht*H>10%1, we may assume that a large fraction will have energies. 
such that 10°'<H'?H<10*!, In this case, from fig. 2, the observed 
ratio, F, will be given by 

0-5< en x 10% <1-6, 
and we may, without much error, put 


FPO A te ic. a ee 
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For hydrogen A is of order 10?6N~-1 cm, so that, if H has the maximum 
value allowed by relation 1, 

gel Ooo Sie ge i8e Grab) mel ea mae) er (6) 

The observed value of this ratio is of order 3x 10-4, if we take the 


y component to contain 0-1°% of the energy of primary cosmic rays 
(Perlow and Kissinger 1951); thus we must put 
EBONY SEES UT tT Rae ieee Bera 8 
If, for example, we suppose N=10°, then we require 7’ 1027 °K, 
which is lower than the accepted value for interstellar gas (104 °K), but 
not improbable if the region were able to radiate thermal energy through 
the presence of even small quantities of dust. In this case we should 
have to assume the order of magnitude of H10~! gauss, and the lower 
energy limit of the electron flux contributing to the radio noise would 
be H>4Mev. Most of the resultant y quantum energies would thus be 
in the range investigated by Perlow and Kissinger (1951). 


Fig. 2 


10° 10 10° 10 ej 10° 
The relation between electron energy and the ratio of its radiation loss 
at radio noise frequencies to that as y quanta. 

The frequency spectrum of the radio noise could only be deduced if 
both the energy spectrum of the electron flux and the variation both 
of the flux and of H throughout the source regions were known. In the 
absence of this knowledge it is only worth while to show that it would 
not be grossly dissimilar from the observed spectrum. Ney: 

If we assume uniform flux and field, take the very rough approximation 
for the spectrum of fig. 1, 

piv)~vi8 for vr, 
p(v)=90 for v>, 
and assume the energy spectrum of the electrons to be 
N(£)~E™, 


then a straightforward integration leads to the noise spectrum 


Herne ) 
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If the electron spectrum is similar to that of primary cosmic rays 
2<a<2-5 and —0-5>(1—a/2)>—0-75, which is compatible with the 
observed spectrum (Ryle 1950). 


§ 3. CONCLUSION 


These expressions must be regarded as very approximate, but they 
serve to show that we should expect considerable variation of source 
strength with temperature and a smaller variation with density. The 
observation that the few discrete radio noise sources at present identified 
with extra-galactic nebuli (Hanbury-Brown and Hazard 1952) show 
some tendency to be associated with those of low optical brightness may 
thus favour the assumption of an inverse relation between magnetic 
field and temperature, such as we have assumed. 

It is clear that even if the magnetic fields produced by the random 
motion of matter are much less than would be required by equipartition 
of energy, plausible conditions can still be postulated in which the 
production of cosmic rays by an electromagnetic process would lead 
naturally to the observed radio noise flux. If such a mechanism is 
indeed an important source of radio noise, the y component of primary 
cosmic rays should show the same marked anisotropy as that of the 
galactic radio noise. More experimental evidence on the magnitude, 
spectrum and direction of this component would be very desirable. 

The presence of nuclei heavier than protons in the primary cosmic- 
ray flux could not be explained by the escape mechanism discussed 
above. If such nuclei were accelerated in regions similar to those 
proposed they would have to escape as charged particles. It would 
probably be an over-simplification, however, to assume that all cosmic 
rays are produced in similar sources. 
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ABSTRACT 


The division of energy between the secondary particles in the direct 
production of electron pairs by fast electrons is discussed, and experimental 
results are compared with the Bhabha theory. The angular separations 
of the tracks of the secondary particles are considered. An event attri- 
buted to pair production by a fast meson is described. 


$1. INTRODUCTION 


A TRIDENT has been defined (Powell 1949) as an event, observed in an 
electron sensitive photographic emulsion, in which an electron pair is 
created on the path of a relativistic charged particle, presumably by its 
interaction with an atomic nucleus. Observations of these events have 
been made in a number of laboratories (Occhialini 1949 ; Bradt, Kaplon 
and Peters 1950). It has been shown (Hooper, King and Morrish 
1951 a) that the great majority of tridents are initiated by fast electrons 
and that they represent a distinct process and are not due to the chance 
conversion of a bremsstrahlung y-ray close to the track of the parent 
particle. In the majority of reported measurements it has been proved 
that, within the limits of experimental error, a balance of energy exists 
between the primary and the secondary particles. 

Bhabha (1935) has considered the theory of the process, and has shown 
that a close relationship exists between pair production by a fast (6~1) 
charged particle and pair creation by a photon. A Fourier analysis may 
be made of the field of the incident particle, so that the components are 
considered as aspectrum of photons. v. Weizsaicker (1934) has shown that 
the form of this spectrum is given approximately by I/hv. Bhabha 
considers the materialization of these photons in the Coulomb field of a 
stationary particle, and thus evaluates the trident cross section by the 
use of the theory of Bethe and Heitler (1934). This treatment leads to an 
expression for the differential cross section of the form 


HOD LAA ai Bat Bid Gadh gem vs if ex. (1) 


where Z,¢ and Z,¢ are the charges of the moving and stationary particles 
and E_, E,, are respectively the energies of the created electron and 
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positron. B represents the relativistic factor [1 —(v?/c?)]-1? for the incident 
particle. Equation (1) holds for values of H_, E.. within the limits : 


mc? <H_, H,<primary energy. 


Integration of (1) with respect to H_, H, leads to an expression which, 
provided the condition 
| B<2-187.Z-'8, = 9 on 


is satisfied, shows that the cross section varies approximately as the cube 
of log,B. For values of B exceeding those given by (2) the effects of 
screening by the atomic electrons must be considered, and Bhabha 
(1935) taking these into account, derives the more complicated form given 
by eqn. (47) in his paper. This, however, leads to a variation not 
markedly different from that given by the simple cubed logarithm relation. 

It has already been shown (Hooper et al. 1951 a) that the experimental 
_ results are compatible with this variation of cross section, but it has not 
been possible to distinguish between various more detailed theories 
(e.g. Racah 1937 ; Ravenhall 1950). 


§2. ENERGY DrvisIon In TRIDENTS 


It is notable that in more than 90° of the tridents observed, one of the 
secondary particles carries away more than half the primary energy ; in 
more than 50°, of the tridents over two thirds of the primary energy is 
carried away by one of the secondaries. It was suggested that this. 
unequal energy division might reduce the otherwise large effect of the 
trident process in the early development of electromagnetic cascades. 
(Hooper et al. 1951 b). 

Integration of (1) over all energies H.,, of the positron leads to an ex- 
pression which shows that the probability of the electron having an energy 
E_, independent of the value of H,, is approximately proportional to 
1/H_. Since (1) is symmetric in H_ and E,, a corresponding relation holds 
for H,, and therefore a distribution in energy of the created electrons 
would, for a particular primary energy, also be of hyperbolic form. 

An important feature of these experiments is our inability to distinguish 
the created electrons with certainty among the secondary particles. Even 
if the charge carried by each particle was measured with a powerful 
magnetic field (Dilworth, Goldsack, Goldschmidt-Clermont, and Levy 
1950), it would be possible to do no more than identify one of the three 
secondaries as a created particle. 

The distribution of the energies of the created electrons, predicted by 
Bhabha, suggests that in the majority of events the two created particles 
are those of lowest energy. In this work we have designated the secondary 
particle energies #,, H,, and H,, in order of decreasing magnitude, so that 
E, and #, are taken, by convention, to represent the materialization of a 
‘ virtual’ photon. 
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Our measurements have been made on 56 tridents which occurred in 
plates exposed to the cosmic radiation ; thus the energy of the primary 
particles may lie between 60 Mev and 100 Bev. Nevertheless, a distribution 
of the values of E/E, where Hy is the primary energy, should be hyperbolic, 
of the form Hy/H—1. A histogram on which the values of both £/Ey 
and #;/Hy are plotted is shown in fig. 1. It may be seen that for the lower 
values of #, ;/Hy the histogram is in satisfactory agreement with the 
hyperbolic curve, which is a normalized plot of the function Zy/E—1. It 
is clear, however, that on account of our convention in designating the 
created particles, the histogram cannot extend beyond HF, 4/Hy=0-5. 


Fig. 1 


No. of Particles. 


Wi 


Hy, 3/Hy 
The distribution in values of both #,/Hy and H,/Ey. The smooth curve 
represents the function Hy/H—1 drawn to enclose an area equal to that of 
the histogram. 


According to the theory, the secondary particle of highest energy is, in a 
small proportion of the tridents, one of the created electrons. If we were, 
in fact, able to identify the created pair particles among the secondaries 
we might reasonably expect better agreement of the histogram with the 
theoretical curve. 

In fig. 2 we have plotted the distribution of the values of (H,+E3)/Ey 
for all the measured tridents. A recent detailed calculation by Ott (1952) 
on the basis of the theory of Bhabha predicts that such a distribution has a 
maximum when the energy of the created pair is about 2 Mev, whatever the 
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energy of the primary electron. For higher energies of the created pair, 
the distribution is again approximately hyperbolic, as would be expected. 
As our primary energies are all greater than 60 Mev, the maximum in the 
curve which Ott predicts does not suffice to explain the departure of the 
histogram from the hyperbolic form. We are unable to attribute this 
result to the use of our convention, although it may be due to poor statistics 
or a failure to observe those tridents in which the created pair is of very low 


energy.* 


Fig. 2 
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The distribution in the values of the quantity (H,+#3)/Hy for each of 56 
tridents. 


Our results in fig. 1 indicate that in the production of tridents by 
electrons one of the created pair particles will be, in about 10% of the 
events, the secondary particle of greatest energy. This figure agrees 
roughly with that given by the Bhabha theory. A more accurate check 
on the theory will be possible when the new high energy electron 
accelerators come into use. 
ee eee 

* Note added in proof :—Ott has i : 
the experimental fetes chown te enau sear pre ewe ek 
more detailed theoretical treatment. He states that this comparison aes 


suggests that a considerable number of trid i 
secondaries have not been observed, ep malay si og Seal 
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§3. ANGULAR SEPARATION IN TRIDENTS 


Bhabha shows further that, with electron pairs directly created by an 
electron, the angular separations of the secondary particles with respect 
to the primary direction are of the order m ,c?/H_, and m,c?/E., respectively. 
On the other hand, in electron pairs created by heavier particles, the 
angular separation is of the order of 1/B, depending on the energy of 
the primary particle. The cross section for pair creation by mesons or 
protons is small compared with that by electrons. 

The small angular separations observed for secondary particles in tridents 
confirm the earlier conclusion (Hooper e¢ al. 1951 a) that the great majority 
of these events are caused by electrons. There are only four events in 
which the angular separation of one of the secondary particles is of the 
order of 5°, but in these examples we have found that the energy of the 
particle concerned was less than 10 Mev, not very large in relation to m,c?. 
A similar event, described by Occhialini (1949) may be explained in this 
way. We therefore consider that our observations lend strong swpport 
to the Bhabha theory. 


Fig. 3 
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A photomosaic of an event in which a singly charged particle (track pp) appears 
to produce an electron pair (ce) at T. It is probable that the incident 
particle was a meson. 


In conclusion we describe in detail a trident which is remarkable in 
several respects. This event is shown in fig. 3. The grain densities in all 
the tracks have the minimum value; scattering measurements on the 
track of the incident particle p indicate a primary energy of about 1200 Mev 
and there is no measurable angular deviation at Tin pp. The two tracks 
marked ¢ are each of energy about 50 Mev, so that it is reasonable to assume 
that ee represent a created electron pair. 

The notable features of the event are the wide angle and the symmetry 
in both angle and energy of the secondary particles. We note that for 
this example the condition 


m c?<H_, H,<primary energy 
is valid, irrespective of the mass assumed for the incident particle. 
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The magnitude of the separation angles (8° each) is consistent with the 
view that the primary particle is a 7 or « meson, for which the value of B 
is about 8. 

It is significant that in the observations on the 5 mm track of the primary 
particle there is no evidence, such as radiative energy loss, which would 
conflict with this interpretation. 
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ABSTRACT 


The Thomas—Fermi statistical model of the atom leads to the formula 
20:92 Z’? ev for the total binding energy of an atom with atomic 
number Z, but this formula gives values which are too high by roughly 
20%. The actual binding energies increase quite smoothly with 
increasing Z, which suggests the existence of a more appropriate formula. 

The reasons for the discrepancy are explained and a new theoretical 
formula derived. This is compared with spectroscopic and self-consistent 
field data and is found to agree within 1 or 2%. Some uncertainty 
remains about the relativistic correction for heavy atoms. 


§1. [INTRODUCTION 


Tue Thomas—Fermi statistical model of the atom leads to the formula 
20:92 Z8 ev for the total binding energy of an atom with atomic 
number Z. It is known that the values given by this formula are some 
20 or 30% higher than the experimental binding energies derived 
from spectroscopy. The actual binding energies of light atoms are 
found to increase very smoothly almost exactly as Z’?, but with a 
different constant (Young 1929, Allard 1948) : 


Vie oedsrZ i wey el Peete wes te. (1) 


A different empirical formula has been proposed by Foldy (1951), which 
is equivalent to 

Weatle- G02 eye te Oe! So (2) 
and which is designed to represent the properties of heavier atoms better. 

The purpose of this paper is to explain the reasons for the discrepancy, 
and to derive a better theoretical formula. This will be shown to be in 
good agreement with the data. 

The failure of the ordinary theoretical formula (due essentially to 
Milne 1927) can be explained as follows. If one seeks to account for the 
total binding energy W of an atom within 1 or 2% by a statistical 
model, it is evidently necessary to make sure that no terms are 
neglected which are larger in order of magnitude than 1/Z times the 
total energy. There are actually two such terms: one arises from 
a boundary effect, and is of the order Z* and the other is due to 
exchange and is of the order 2°”. 


* Communicated by the Author. 
SER. 7, VOL. 43, NO. 343.—AUGUST 1952 3N 
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It will be shown that these terms, which have different signs, combine 
to produce a variation roughly as Z’’* for light atoms (but with a reduced 
coefficient), in agreement with the data of Young and Allard, and a 
variation roughly as Z2/5 for heavier atoms, in agreement with Foldy’s 
data. 

§2. PRELIMINARIES 
2.1. Notation 

All quantities will be expressed in Hartree’s atomic units, unless 
otherwise specified. Following substantially the notation of Schiff (1949) 
and Condon and Shortley (1935), the electrostatic potential at a distance r 
from the nucleus is expressed as a fraction x(x) of the potential due to 
the nucleus itself, where r—bx; 6 is a scale factor proportional to the 
extent of the main part of the charge distribution : ' 


b= (9n?/128Z)8—0-88534Z-N8,  . | (8) 


The total binding energy is denoted by W, and the electrostatic potential 
at the nucleus due to the electrons only is denoted by —v. The symbol v 
will be used in this sense throughout the present paper. 


2.2. Nature of the Statistical Method 

In the statistical method, it is convenient to think of a fictitious or 
purely mathematical problem, in which Z is very large, say a million, 
so that even Z¥8 is large, and the terms of any series in descending 
powers of Z1/8 can be clearly separated. The solutions of a non-relativistic 
wave equation are envisaged, with magnetic effects omitted (i.e. coo, 
a0). The consideration of such an idealized limit helps to clarify the 
procedure, in the present work just as in the original Thomas—Fermi 
calculation, despite the air of unreality which it creates. 


2.3. Connection between W and the Electrostatic Potential of the Electrons 

Of the various ways of calculating the binding energy, perhaps the 
easiest is to relate it to the electrostatic potential —v at the nucleus due 
to the charge distribution of the electrons. The argument consists in 
considering the effect of increasing the nuclear charge by a small 
amount 6Z, This change may be regarded as a perturbing field 
V(r)=—68Z/r. According to standard perturbation theory, the change 
in the energy of the ground state is given by an integral containing the 
perturbing potential and the unperturbed wave functions : 


J x Vy=—s2 [yep 
But this is equal to the potential at the origin due to a charge-density 
—6Z|%|?, and is equal to —v defined earlier. Thus 


dw 
adZ =v. eS ee eee a ene (4) 
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The energy of any electrons which have to be added on the outside of 
the atom to keep it electrically neutral has been here neglected, but the 
correction is O(1) which is smaller than any terms which need be retained 
in the subsequent calculations. 

The usual formula for the binding energy is obtained by noting that 


v=lim (Z/r)(1—y) 
r—>0 


=—(Z/b) . x’(0). Serene, are MO) 
Here 6 is given by (3), and the derivative —y’(0) is Baker’s constant 
er seta 35 las eee Sd ce te) (GO) 
This leads to the usually-quoted formulae 
v=1-794 Z48 eI ies APL CT) 
(Lamb 1951, Mack 1950) and, using (4), 
Wie O60 Ze 
=20;92 2° electron-volts. so. 4 (8) 


The latter formula was first given by Milne (1927), but with a numerical 
inaccuracy which was corrected by subsequent writers, e.g. Baker (1930). 

Those corrections which are larger than 1/Z times the main term (7) 
or (8), as regards order of magnitude, will now be investigated. 


§3. THE Bounpary EFFECT 

The boundary effect, which will be considered next, is perhaps most 
familiar for the case of independent particles confined in a large cubical 
box (volume L*). In the zero-order approximation there are L* .4aP%/h3 
single-particle wave functions with momentum less than P; and this 
~ ean be equated to the number JN of particles when one particle is assigned 
to each state. The next approximation, which is easily obtained, 
corresponds physically to the fact that the density of particles is 
considerably reduced in the region near the boundary within a distance 
of say a quarter of the wavelength of the fastest particle, while at 
interior points the density of particles must be increased by a factor 
1+O(N-1*) so as to compensate. 

In the present application, there are two important differences. First, 
the boundary which concerns us most is a point instead of a surface. 
Although the zero-order (Thomas—Fermi) formulae give an electron- 
density which tends to infinity at the nucleus, it can be shown that the 
necessary boundary corrections reduce the number of electrons in the 
boundary region by a finite amount (finite as Z-> oo), and that the density 
elsewhere is increased only by a factor 1+O(Z~'), which is negligible 
for the purpose of the present calculations. Secondly, the field is 
singular, so that a finite error in the number of electrons in the boundary 
region (radius O(1/Z), or roughly the K-shell region) causes an error of 
the order of Z? in the binding energy, that is to say a fraction 1/Z' of 
the total binding energy, which cannot be ignored in practice. 


3N2 
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To calculate this correction quantitatively, the true effect of the 
K and L shells will be compared with that which is implied by the 
statistical model. It turns out that hydrogen wavefunctions may be 
used because the field in the boundary region is sufficiently nearly a 
Coulomb potential. In the idealized Thomas—Fermi model, the total 
number of electrons with energy less than —Z, within a spherical shell 
of thickness dr, is easily shown to be 


Z (x OHNMARSS 
= |(-— > 3 '< .o $6 @gaetgor ed ae 
Bs (2 7) y” dr. : (9) 
Since a shell of charge Q contributes an amount Q/r to the potential (—v) 


at the nucleus, the electrons with energy less than —H ought to 
contribute an amount 


p [qx zp ett ae oe ee 


to the potential —v at the nucleus. 
Now choose £ so that there are just two electrons beyond this energy 
according to the statistical model. The contribution to v is found to be 


pee ta Gh oy L 


The actual potential due to the two K-electrons on the other hand is 
(minus) 2Z atomic units. Thus the value of —v obtained statistically 
needs to be increased by an amount 


2Z(31/3—1)=0-88Z. 


(In deriving this result y has been replaced by 1, but the percentage 
errors thereby committed are of the order of Z~*/3, which is within the 
accuracy aimed at.) 

If the eight L-electrons are included as well, a similar calculation 
shows that the discrepancy for the first ten electrons is 


2Z(151/32)—0-93Z. 


In the idealized problem in which Z is so enormous that many shells 
of tightly-bound electrons can be taken into account in this way, the 
complete boundary correction to be applied to v is 


—2Z lim {[$n(n+1)(2n+1)]12®—n}=—Z. 


The corresponding correction to be applied to the total binding 
energy W is found from (4) and is 


—3Z’, ©, ALIN S Nua) oe eee 

It is negative because the statistical model allows too great a density of 
electrons very close to the nucleus. 

(It may be remarked that the number of electrons of energy <—EH 


in a Coulomb field is given by the statistical method as well as can 
reasonably be expected ; it is in v that it gives a systematic error.) 
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§4. Toe Excuancre Errect 


The exchange energy, that is to say the contribution to the total 
energy which appears in the form of exchange integrals when we 
construct an approximate wave function out of combinations of 
one-particle wave functions with the correct symmetry properties, tends 
to lower the energy of the ground state. Exchange is disregarded in the 
Thomas—Fermi model, which corresponds to the original Hartree equations 
rather than to the Hartree-Fock equations. Its effect in a statistical 
model has been discussed by Dirac (1930). 

It will not be necessary to report our calculations in detail, as a 
textbook has since appeared (Corson 1951) which covers most of the 
essential points. It should, however, be added that it is not necessary 
to abandon the Thomas—Fermi distribution and calculate a new 
distribution allowing for exchange (this in any case leads to serious 
difficulties in the peripheral region of the atom). The reason is that the 
effect of exchange is to modify the density distribution by a fraction 
O(Z-*'8), except at the extreme outside of the atom which does not affect 
the total energy appreciably ; and the energy may be calculated from 
the unperturbed density-distribution with only a second-order error 
provided that an appropriate method is used, as is discussed in the usual 
accounts of variational methods. 

The result (cf. Corson 1951, eqn. (9.48), p. 169) is a positive contribution 
to the binding energy W, of magnitude 


2 ro 
=) Na et AE deo Liv eh ce) 4(12) 
0 


7 


where 6 is given by (3), and the integral has the value 
{ Rie 0:6) Gilg on be ied) (13) 
0 


(It was calculated from the table of Bush and Caldwell (1931), which is 
only accurate to $%.)* 


§5. RESULTS 


Collecting the results, it is found from (8), (11), (12) and (13) that the 
total binding energy of an atom of atomic number Z, in the non- 
relativistic approximation, is 


W=0-7687Z78 —1Z7°+ 0:221Z5810(Z48) . . . (14a) 
in Hartree units, or in electron volts : 
20-92Z7/8— 13-60Z2+ 6-01Z5/24+ O(Z4/8) ev. . . . (146) 


For the potential at the nucleus due to the electrostatic field of the 
electrons, we find : 
y= 1-7937Z48—Z-+0-368Z7/24O(Z28). . . 2. . (15) 
Peter cher ahs Of. he Ss OS ee 
* Note added in proof : The constant (13) is equal to 0°61543. 
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$6. DISCUSSION 


6.1. Comparison with Observation 
For the first twenty elements or so, much accurate data on ionization 
energies has been accumulated by spectroscopists, even for the high 
stages of ionization which concern us here (Moore 1949). The last two 
electrons (K shell) constitute a special case, and the energy required to 
remove them has been obtained from the very accurate formula of 
Hylleraas (1930), which being non-relativistic is more directly comparable 
with the present results than are the actual observations (which indeed 
do not extend beyond fluorine); the relativistic correction (see e.g. 
Bethe (1933), p. 383) need not be considered at this stage. For the 
L shell, the values reported in the literature extend as far as phosphorus 
(Z=15); and though in some cases they are extrapolations by the 
spectroscopists, they are quite accurate enough for the present purpose. 
With heavier elements, the energy needed to remove the whole L shell 
(ignoring relativity) may be taken from the semi-empirical formula 
(Z—S/?+C. 12 (16) 
Here S has the theoretical value 4+ 69/256—430/2187=4-073 ; and the 
value of C is about —0-2, based on the sum of values from a study of the 
higher ionization potentials in each of the eight iso-electronic sequences 
(cf. Scott 1952). For the M shell, the last-mentioned method has been 


used to supply the ionization potentials of vanadium IX, X and XI; 
the other values used are observational. 


Table 1. Total ‘ Non-relativistic ’ Binding Energy 


K_ shell 
Ls shell 
M shell 
Total 


Cale’d (eqn. 14) 
Cale’d (T-F-M) 


The binding energies W, measured or estimated in this way, which 
should be quite reliable, are compared in table 1 with the Thomas—Fermi— 
Milne formula (8) and with the new formula (14). It is seen that the 
error decreases smoothly from 3% at carbon to 13% at vanadium. 
Since the formula is theoretically only accurate to 1 part in Z, this is 
surprisingly close agreement. 


6.2. Comparison with Self-consistent Field Data 


For the heavier elements, although there are no relevant experimental 
data, it is possible to check the formulae against the results of 
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self-consistent field computations. As the data beyond copper (Z=29) 
do not include the effect of exchange, they have been adjusted by means 
of the semi-empirical formula 


0-369Z?/3— 0-64Z1/8 Sievetce se hoe | (LZ) 


which represents the difference between the Hartree-Fock and the 
original Hartree values satisfactorily in those cases in which both are 
known; the correction amounts to 4 units for mercury (one would not 
expect it to be as large as the full ‘ exchange term ’ for the reasons given 
by Slater 1951). The first coefficient in (17) is the theoretical constant 
used in (15), 

A representative selection of self-consistent field data is given in 
table 2, based on Dickinson’s table, for comparison with the values given 
by (15) and those obtained from the older formula (7) which 
corresponds to the one-term Thomas—Fermi—Milne formula (8). The 
agreement is within 1%, which again is better than we have any right to 
expect from a calculation accurate to 1 part in Z. 


Table 2. Values of v 


Calculated 
(T-F-M) (eqn. 15) 


24-02 18-37 
33°58 26-18 
43-88 34-71 


84-61 69-15 
97-38 80-10 
138-16 115-40 


167-2 140-8 
189-9 160-7 


557 490 
618 545 


6.3. Comparison with Empirical Formulae 

Besides making the foregoing comparisons, which provide a check ' 
of (14) and (15), it is interesting to examine the power laws which have 
been advocated, and which represent the available information fairly 
successfully. From the present point of view, the question is: what 
power law would most nearly coincide with (14) and (15) ? 

Foldy (1951) has plotted the available values of log v against log Z, 
and obtained a power law vx Z” to represent the data between Z=4 
and Z—80, with an exponent n=1-4. The relation (15) is of course not 
a simple power law, but it is natural to inquire what the gradient of a 
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log-log plot of this kind would be when log Z lies in this region. Taking 
the geometric mean of 4 and 80, Z=18, as a representative value, it is 


found that 


d(logv) 4 2 —2/8 
2 = ~ + 0-186Z-18— 0-0383Z-*/ 
d(logZ) 3 a 


—=1:399. 2). ah) oil alee nena 


Thus v varies approximately as the (7/5)th power of Z in this region. 
This empirical law ought in fact to be rather good, except for the heaviest 
elements ; and with them, although it may deviate from the theoretical 
relation (15), both relations will require appreciable relativity corrections. 

The success of the 7/3 power law (1) for light elements arises from the 
fact that the deviation produced by the second term in (14) is partly 
compensated in this range by the effect of the third term, and it seems 
that the neglected fourth term acts in the opposite sense too. Thus, 
over a limited range, the combined effect of these terms is nearly a 
constant fraction of the main term. The coefficient will evidently be 
about 15 or 16 ev, e.g. if we put Z=8 we find 


W/Z7/?=20-92—4 x 13-60+4X6-01=15-62ev . . . (19) 


to be compared with 15:68 (Young 1929), 15:73 (Allard 1948) or 
15-6 (Condon and Shortley 1935) for the elements He to F. But for the 
very lightest elements the statistical method is so inappropriate that any 
success it has must be ascribed mainly to chance. And for the medium 
and heavy elements (1) is certainly wrong. 


§7. RELATIVITY AND MaGnetic EFFECTS 


Relativity effects of all kinds have been disregarded so far. Though 
this simplification has no serious consequences for Z<30, these effects 
are quite important for heavy elements. It would be a difficult task to 
calculate them accurately. A straightforward extension of Thomas’ 
statistical method (Vallarta and Rosen 1932) is inapplicable to our 
present problem, because most of the correction originates in the region 
close to the nucleus where the statistical method is vitiated by the 
‘boundary effect’, and, in fact, such,methods would give an infinite 
binding energy. Moreover, the interaction between the electrons is 
not wholly electrostatic. 

The relativistic correction to the energy required to strip off the last 
two electrons is, however, accurately known, and the correction for 
stripping off all the others is only of the same order of magnitude. We 
have made rough numerical estimates of the necessary correction to W in 
various ways, and find that our results can be conveniently represented 
by the approximate expression 


Wr- W wr=4 x 10-8 Z9/2 . . e e e ° (20) 
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in atomic units. It amounts to a fraction of a unit as far as phosphorus, 
and increases from about 15 at Cu to two or three thousand at U. The 
error in the total binding energy resulting from neglecting relativity is 
roughly (Z/30)?%. 

§8. CONCLUSION 


The total binding energy of an atom neglecting relativity effects is 
given by (14), with an error which is in principle of the order of magnitude 
of 1 part in Z, and which seems in fact to be about 1 part in 4Z. An 
estimate of the relativity correction is given by (20). The failure of the 
currently-quoted formula is due partly to the shortcomings of the 
statistical model in the region nearest the nucleus, and partly to the 
effect of exchange. 
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SUMMARY 


A quantum mechanical calculation of the lattice constant, 
compressibility, cohesive energy and work function of metallic beryllium 
is carried out, using an extension of the cellular method of Wigner and 
Seitz. Wave functions and energies corresponding to the lowest 
electronic state in the metal are obtained by numerical integration of 
the Schrédinger equation, and the ion-core field used is of the Hartree 
type. The Fermi energy is evaluated by determining the effective 
electronic mass according to the method devised by Bardeen. Free- 
electron values are used for the exchange and correlation energies ; 
it is pointed out that the corrections to these two terms may tend to cancel 
each other. The energy of the valence electrons in the free atom is 
obtained from the calculations of Hartree, together with the correlation 
energy. With the exception of the work function, the final results agree 
very satisfactorily with experiment and compare favourably with the 
values obtained by more elaborate methods. 


$1. INTRODUCTION 


THE quantum theory of metals, as formulated in the cellular approxima 
tion of Wigner and Seitz (1933, 1934), has proved eminently successful 
in the calculation of the cohesive energies of the monovalent metals 
but as yet has not been applied extensively to metals of higher valency. 
The first comparable calculation for a polyvalent metal was carried 
out by Herring and Hill (1940) for beryllium ; these authors determined 
directly the energies at a number of points in k-space and thus obtained 
the density of states and the mean Fermi energy. More recently, a 
simpler method, related more directly to the Wigner—Seitz treatment, 
has been applied by Raimes (1950) to magnesium. 

The present paper contains the results of a new calculation for 
beryllium, similar to that of Raimes, undertaken with the object of 
ascertaining whether a simpler treatment than that of Herring and Hill 
would be satisfactory in this case, and whether reasonable results could 
be obtained without taking explicit account of the zone boundaries. 
Furthermore it was hoped that additional evidence could be obtained 


* Communicated by Professor W. Sucksmith, F.R.S. 
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regarding Raimes’ estimate of the correlation energy in the free atom, 
which consists of taking a proportion of the total energy of the valence 
electrons or of the Coulomb interaction energy. 

A brief note summarizing the results of the present calculation has 
already been published (Donovan 1951). 


§2. CALCULATION OF WaAvE FUNCTIONS 


Metallic beryllium has a close-packed hexagonal structure, and we 
shall assume that it is a reasonable approximation to replace the atomic 
polyhedron for this lattice by a sphere of equal volume. In constructing 
a self-consistent field in the metal we neglect exchange interaction 
between the valence and core electrons ; this will be neglected also when 
considering the free atom. Moreover, in the derivation of the total 
energy of the metal, exchange and van der Waals interactions between 
the ion-cores are neglected.* 

The self-consistent field we shall use is based on the fact that the 
probability of two electrons with parallel spins being in the same atomic 
cell is very small. Thus for a given electron the potential field may be 
expressed as the sum of the ion-core field and that due to the charge 
distribution, within the cell, of all the electrons with opposite spin 
(cf. Raimes, loc. cit., eqn. (3)). The ion-core field was taken from the 
Hartree calculation for Bet++ (Hartree and Hartree 1935 a), assuming 
that the ion-core remains unchanged in the metal. 

The wave function y%, corresponding to the lowest electronic state, 
and the eigenvalue «, of the appropriate Schrédinger equation, may be 
obtained by the iterative process of numerical integration described by 
Raimes (loc. cit.). This procedure was carried out for eight values of 
the radius 7, of the atomic sphere and the resulting curve for ¢) is shown 
in fig. 3. The general form of the wave functions may be seen from 
fig. 1, which represents y% for r,=2-30 Bohr units. This wave function 
is normalized according to the equation 


Ts 
[ “12.2 dr=1. 
“0 


$3. Tora, ENERGY IN THE METAL 


On the basis of the assumptions mentioned in § 2, the total energy of 
the valence electrons in the metal may be calculated according to the 
one-electron approximation. As a consequence of the particular 
potential used the resulting expression does not include a separate term 
for the Coulomb energy (cf. Raimes, loc. cit.) and the total energy, 
expressed in rydbergs per electron, is given by 

a(r,) 1154 0-576 
ake Sip ye NC ye cr ee eee CL 
eke ae aes r, 0-797, +51 (1) 
i 

* The validity of this assumption has been verified by Herring and Hill 
(loc. cat.). 
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In this expression « is the effective mass parameter, i.e. the ratio of the 
actual electronic mass to the effective mass near the centre of the 
Brillouin zone. The second, third and fourth terms in (1) represent 
respectively the Fermi, exchange and correlation energies. 

In writing the Fermi energy as the product of « and the free-electron 
value we are tacitly assuming that all the electrons have the same 
effective mass and we are thus ignoring the effect of the zone boundaries. 
Using the method devised by Bardeen (1938) the value of « was 
determined for every value of r, and the resulting curve is shown in 
fig. 2. The variation with r, closely resembles that for lithium and 
Bardeen’s values for this case are plotted in fig. 2 for comparison. 


Fig. 1 
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Wave function for lowest electronic state, corresponding to r =2-30 Bohr units. 


Herring and Hill (Joc. cit.) calculated values of « for three values of r,, 
as indicated in fig. 2, and it is interesting that their value of 0-616 for 
,=2-37 is reproduced exactly by the present work. However a curve 
drawn through the three points would lie slightly above our fig. 2 for 
r,<2-37 and considerably below for r,> 2-37. In this connection it should 
be pointed out that in the calculations of Herring and Hill the same 
potential (that for r,—2-37) was used in all three cases whereas in the 
present treatment each integration for « was carried out with the 
potential function appropriate to the particular value of t 

As indicated in (1) free-electron formulae are used for the exchange 
and correlation energies, discussed respectively by Bloch (1929) and 
Wigner (1934). Tentative suggestions for modifying the free-electron 
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formulae have been put forward (Herring and Hill, loc. cit., Herring 1951) 
but, on applying these to the present calculation, it is found that the 
corrections to the exchange and correlation energies cancel each other in 
such a way as to leave the sum of these terms virtually the same as for 
free electrons. 

The values of H obtained from (1) are plotted in fig. 3. If the 
free-electron expression be used for the Fermi energy the values of the 
total energy do not lead to results with quantitative significance. This 


Fig. 2 


6 20 24 23 a2 36 
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Effective mass parameter for beryllium and lithium. The circles represent the 


values calculated by Herring and Hill. The broken vertical lines indicate 
the observed lattice constants. 


behaviour is in contrast to that for magnesium, for which Raimes 
obtained better agreement with experiment by using the free-electron 
value for the Fermi energy. 


§4. ToraL ENERGY IN THE FREE ATOM 


For the free neutral atom the appropriate data may be obtained from 
the calculations of Hartree and Hartree (1935a, b). Since we are 
neglecting exchange interaction between the valence and core electrons 
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both in the metal and in the free atom, we require firstly the energy of 
the valence electrons according to the Hartree approximation. This is 
—0-940 rydbergs per electron (line a in fig. 3). 

Since the two valence electrons have opposite spins there will be no 
exchange energy, but some cognisance must be taken of the correlation 
energy. This is a measure of the error involved in the use of one-electron 
wave functions and should therefore be given by the difference between 
the experimental and the Fock values of the total energy. For the 
2s electrons the correlation energy can be estimated only in an indirect 
manner and the position is unsatisfactory in view of the discrepancies 
between the various estimates. 


Fig. 3 
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Kigenvalue of lowest electronic state and total energy in the metal. The lines 
a and 6 represent respectively the total energy in the free atom with and 
without correlation energy. 


Raimes (loc. cit.) has suggested that the correlation energy might be 
estimated by taking a definite proportion of the observed total energy 
or of the Coulomb energy for the relevant electrons. He obtained the 
numerical factors from the helium atom, for which the total correlation 
energy is known to be —0-076 rydbergs. On this basis the correlation 
energy 1s estimated by multiplying the observed total energy by 0-013 or 


of Some Properties of Metallic Beryllium 873 


the Coulomb energy by —0-037. For beryllium both methods give the 
same result, namely, —0-013 rydbergs per electron. If this value be 
accepted the total energy in the free atom is —0-953 rydbergs per electron 
and this is indicated by line 6 in fig. 3. 

However, an estimate of —0-05 rydbergs per electron for the 
correlation energy has been given by Seitz (1940, p. 248) and it is 
difficult to reconcile Raimes’ suggestion with the available data for 
beryllium. In particular, it leads to a much larger correlation for the 
1s electrons than for the 2s electrons, which is extremely unlikely in 
view of the fact that the total correlation energy in Be is more than 
double that in Bet+. 


§ 5. Resutts anp Discussion 


_ The equilibrium lattice constant, compressibility and cohesive energy 
may be determined by the usual methods and the results are set out in 
table 1, together with the experimental values and those obtained by 
Herring and Hill. In the latter case extreme values are quoted where 
various methods were employed in the calculations. The experimental 
value for the lattice constant may be deduced from the measurements 
of Gordon (1949) and the experimental compressibility has been given 
by Bridgman (1932). 


Table 1 
| 
Lattice Compressi- | Cohesive Work 
Constant bility Energy | Function 

(Bohr units) | (cm?/kg x 107) |(K cal/mole)| (volts) 

Present 
Calculation 2-34 8-0 84 —0:3 
Herring and Hill eco Toles < Alihog 53 —1-9 
Experimental 2-36 7:8 75 3°9 


(0° c) (30° c) 


The cohesive energy is given by the difference in energy between 
line 6 and the minimum of # in fig. 3 and the rather high value obtained 
tends to strengthen the view that the correlation energy has been 
under-estimated by Raimes’ method. If the correlation energy is taken 
to be —0-05 rydbergs per electron, as mentioned in § 4, the binding 
energy is 61K cal/mole. Herring and Hill used the experimental value 
for the energy of the valence electrons which, as they admitted, is too 
large and leads to a value of the cohesive energy which is much too low. 

Finally, a rough value for the work function may be obtained using 
the formulation of Herring and Nichols (1949, p. 230). The value given 
in table 1, which excludes the contribution due to the surface double 
layer, was derived from the value of # at the calculated lattice constant. 
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Comparison with the observed value (Mann and DuBridge 1937) shows 
that the disagreement, though smaller than that found by Herring and 
Hill, is still too large to be explained by the double layer term. 

With the exception of the work function, the calculated values obtained 
in the present work are in good agreement with experiment and are 
certainly no less satisfactory than those of Herring and Hill. However, 
it must be emphasized that in the present application, it is necessary to 
consider variations only in the size of the atomic polyhedron ; the form 
of the cell remains unaltered.. In view of the neglect of the zone 
boundaries the present method would be of little value in problems 
involving a distortion of the atomic polyhedron. 


The author wishes to express grateful thanks to Professor A. G. Walker 
and to Mr. Ge H. Jowett for placing at his disposal a Marchant electric. 
calculating machine, on which the numerical work was carried out. 
The author is also indebted to Dr. 8. Raimes, of Imperial College, 
London, for helpful conversations. 
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ABSTRACT 


Soft electromagnetic radiations from 22®Pu and 2°3U have been investi- 
gated with a proportional counter filled with krypton or xenon. 

*89Pu emits the L x-radiation of uranium with an intensity of 
(4+1)x10 L x-rays per « particle. The L x-radiation arises as a result 
of the internal conversion of y rays. Residual unconverted gamma 
rays with energies of 52-0+0-3 kev and 38-5-L0-4 kev are emitted. The 
intensities of these gamma rays are 7x 10-5 and 210-5 quanta per « 
particle respectively. K x-rays are emitted with an intensity of about 
2x 10-° quanta per « particle. 

A gamma radiation of energy equal to 45-0+0-3 kev is ascribed to 
pb or.-*°Pu. 

283 emits the L x-radiation of thorium with an intensity of (4-1) x 10-? 
L rays per « particle. Gamma radiations of energies 42-8--0-3 kev and 
56-1+0-4 kev are present with intensities of 510-4 and 10-* quanta 
per « particle respectively. 


§1. INTRODUCTION 


FINE structure is present in the energy spectrum of « particles from 
many of the natural « particle emitters. The fine structure groups arise 
from transitions to excited states of the daughter nucleus and so their 
x particles are of lower energy than those which arise from the transition 
to the ground state. The detection of fine structure in the « particle 
energy spectrum itself is often difficult due to the small energy differences 
involved and to the usually weak intensity of the fine structure components. 
In these cases the radiations emitted from the excited states of the 
daughter nucleus following the emission of an « particle are often more 
readily detectable and provide much more accurate values of the energy 
levels of the residual nucleus than does a direct measurement of the 
x particle spectrum. — 


* Some of the results obtained with plutonium have already been briefly 
reported (West and Dawson 1951). 
+ Communicated by the Authors. 
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The proportional counter is a very convenient instrument for investi- 
gating y radiations in the energy region up to about 100 kev. Filled 
with krypton or xenon, such counters have a high efficiency for detecting 
y radiations ; moreover the great majority of the y radiation detected 
is absorbed in the gas of the counter by the photo-electric process, so that 
the size of pulse produced in the counter is directly proportional to the 
total energy of the y radiation. (Hanna, Kirkwood and Pontecorvo 
1949.) 

Soft radiations from 72°Pu were first investigated by Ghiorso (1944) 
by an absorption method. Rosenblum, Valadares and Goldschmidt 
(1950) recently reported the presence of fine structure in the « particle 
energy spectrum of *8*Pu. They observed a low energy group of alpha 
particles differing in energy by about 50 kev from the main group. 
The present work was undertaken to examine the radiations from the 
excited state of ?°°U which results from the emission of a low energy 
group of « particles by *°*Pu. The method is of general application and 
results obtained with 7°°U are also described. 


§2. APPARATUS AND METHOD 


The proportional counter was made of soft glass 1 mm thick with an 
aquadag cathode on the outside (Maze 1946). Its diameter was 5 cm 
and the active length was 34cm. It was fitted with a tungsten central 
wire of diameter 0-005”. Half way along the counter two thin mica 
windows (0-0038 cm thickness) were incorporated at opposite ends of a 
diameter to admit very soft radiations. The counter was filled to about 
55 cm of mercury pressure with a mixture of krypton (or xenon) and 
10% methane. Pulses from the counter were amplified by a head 
amplifier and linear amplifier, the output of which was fed into a 30: 
channel pulse analyser whose channels were | volt wide. , 

Energy calibrations of the counter were carried out with characteristic 
K x-radiation from K capture bodies (Rothwell and West 1950). Calibra- 
tion runs were taken immediately before and after the measurement 
with the radiation under investigation to eliminate the effects of any 
slight drifts in the apparatus. Any lack of linearity in the amplifier 
and the associated electronic equipment was eliminated by calibrating the 
apparatus with artificial pulses from a pulse generator. The size of pulse 
fed into the apparatus could be varied by accurately known amounts 
and a curve of output pulse size versus input pulse size was plotted. All 
pulse sizes are then referred by means of this calibration graph to the 
output of the pulse generator which was known to be linear. 

Strong « particle sources excite the characteristic x-rays of materials 
which they strike (Curie and Joliot 1931). For this reason sources were — 
allowed in contact only with light elements whose characteristic X-rays 
are very soft. Where a thin source was required it was deposited on 
aluminium ; stronger sources consisted of solid material contained in 
thin glass tubes. 
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§ 3. ErrictIency or THE CounTER ror y RapIATIoNs 


In the absence of any escape of photo-electrons from the gas into the 
walls of the counter, the efficiency of the counter can be calculated from 
the geometry and the photo-electric absorption coefficient of the counter 
gas for the y radiation. The loss of photo-electrons due to wall effect 
was investigated by varying the pressure of krypton in the counter. 
Uncollimated sources of ¥Gd (41-4 kev), 244Am (60 kev) and 4°Tm 
(85 kev) were counted in a standard position in contact with the outside 
wall of the counter. A typical pulse size distribution from 1Gd is 
shown in fig. 1. The ‘escape’ peaks having an energy 12-7 kev below 


Fig. 1 
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that of the peaks corresponding to the full energy of the radiations is 
characteristic of krypton (West and Rothwell 1950). The escape peak 
occurs as a result of the escape of the K x-radiation of krypton (12-7 kev) 
from the counter following photo-electric absorption in the K shell. 

The counting rate at a given pressure of gas was mea ured for each 
of the sources by evaluating the areas under the ‘ escape ’ and full energy 
peaks in the pulse size distributions. The mean free path for photo- 
electric absorption is large compared with the counter dimensions at the 
energies and pressures in question, so the number of quanta absorbed is 


402 
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proportional to the pressure of krypton in the counter. The counting 
rate per unit pressure is therefore a measure of the efficiency of detecting 
an absorbed quantum. The pressure of krypton could not be made large 
enough to obtain the saturation values for each radiation corresponding 
to no escape of photo-electrons. Instead the counter was placed in a 
magnetic field of 6000 gauss parallel to its axis (West and Rothwell 
1950) and in these conditions it was assumed that the amount of escape 
was negligible. The counting rate per unit pressure of krypton is 
normalized to unity in the presence of the magnetic field. The curve of 
efficiency of detection of an absorbed quantum versus pressure of krypton 
is Shown in fig. 2. 


Fig. 2 
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It was noticed that even when as much as 50°% escape occurred there 
was no appreciable broadening of the peaks in the distribution. The 
photo-electron follows a tortuous path and the wall effect sets in when 
any part of the track grazes the counter wall and not when the end of its 
track just reaches the wall as is the case with particles following a 
rectilinear path. With soft electrons, therefore, the onset of wall effect 
is characterized by an immediate large diminution of pulse size which 
removes this pulse completely away from the region of the peak. The 
fact that a narrow peak is obtained is therefore no indication that the 
wall effect is small. 
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The amount of escape in a krypton filled counter at these energies is 
mainly a property of the energy of the photo-electron removed from the 
K shell. It is possible, therefore, to make an approximate calculation 
of the escape in other gases for radiations which give the same energies 
of photo-electron as those measured in krypton. It is assumed that the 
effective range of the photo-electron is inversely proportional to the 
number of electrons per c.c. in the gas. This assumption will result 
in an underestimation of the efficiency for gases of atomic number greater 
than krypton and an overestimation for gases of smaller atomic number. 
The derived curves for argon and xenon are shown in fig. 3. 


Fig. 3 
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Approximate curves (derived from measurements in krypton) showing efficiency 
of detection of an absorbed y ray in a 5 cm diameter counter filled with 
argon or xenon. 


§4, MEASUREMENTS WITH 7?9Pu 


The L X-radvation : 

A source of 66 micrograms of 789Pu was deposited on aluminium and 
covered with a thin film of ‘ Zapon’ lacquer. It was placed close to one 
of the windows of a counter which contained krypton and the pulse 
distribution measured. The distribution obtained is shown in fig. 4. 
Tt consists of three peaks which are characteristic of atomic L X-radiation. 
The peaks correspond to the «, 6 and y groups of L x-ray lines. The 
energies of the groups were measured from calibrations of the counter 
with 20-2 kev x-rays from a !°%Pd K-capture source (17 days half-life). 
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The escape peaks from the L x-radiation shown in fig. 4 are shown in 
fig. 5. There is no escape peak from the L, radiation which must therefore 
have an energy less than 14:3 kev (the K binding energy of krypton). 
The energy of the escape peaks was measured by calibrating the counter 
with 8-07 kev x-rays from ®Zn (250 days). The energies of the Lg and 
L, radiations could then be obtained by adding 12-7 kev in each case. 
The percentage accuracy of the determination is thereby improved. 
Several independent measurements were made; the results are given 
in table 1 which includes the energies of the L radiations of uranium and. 
plutonium for comparison. 


Fig. 4 


2000}- 


Counts per channel 
ji 


J/OOO}- 


Pulse size (volts) 
L x-radiation from 239Py, 


Table 1* 

L, Ls Ly, 
Experimental values (kev) 13-2+40-3 17-0-++0-2 20-4+0-2 
Energies of uranium L x-rays (kev) 13-6 16-9 20-2 
Energies of plutonium L x-rays (kev) 14-2 17-9 21-4 


*'The L x-ray spectrum is complex and consists of many individual lines 
which are not resolved by the proportional counter technique. The energy 
values quoted are the mean energies of the more prominent lines within each 
group. The values for plutonium are obtained by extrapolation. 
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It is evident from table 1 that the radiation is the L x-radiation of 
uranium, An estimate of its intensity was obtained by measuring the 
counting rate in each group of L rays when the counter was exposed to a 
collimated beam of the radiation. The efficiency of the counter was 
estimated from the geometry and the known absorption coefficients of 
krypton for the individual radiations. The measured intensity was 
4x 10 L rays per « particle with an estimated uncertainty of +1 10-2. 


The L,, Lg and L, radiations are present in the ratios 1:0: 1:5: 0-3 
respectively. 


Fig. 5 
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§ 5. Sorr Gamma RavbIaTIONS 


Very much stronger sources of ***Pu were needed for studying the 
soft gamma radiations. Sources several milligrams in weight were used. 
Each source was placed in contact with the glass wall of the counter 
instead of over the window to reduce the intensity of L x-radiation 
detected. Measurements were first carried out using a counter containing 
krypton. The pulse size distribution is shown in fig. 6. The energies 
of the radiations were determined after the counter had been calibrated 
with characteristic X-rays of energy 41-4 kev from 1°Gd (half-life 155 
days). The interpretation of the pulse size distribution shown in fig. 6 
is complicated by the presence of the escape peaks. The peak A at 
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52:3-10-4 kev energy represents the total energy of one radiation since 
there is no comparable peak at 52-3-+-12-7=65 key. The escape peak 
from this radiation should give rise to a peak of about equal intensity 
at 39-6 kev. The peak B at 38-4+0-4 kev energy is too large to be 
entirely due to the escape peak of the 52-3 kev radiation. A radiation 
of energy close to 38:4 kev must be present as well. The energy of this 
second radiation, however, cannot be determined accurately due to its 
accidental coincidence with the escape peak. The weaker radiations 
consist of peak C at 59-2+0-5 kev, whose escape peak lies at 46-5 kev 
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It is not possible from these measurements alone to decide whether the 
peak at 44-9+0-4 kev (peak D) can be accounted for entirely as an 
escape peak, in which case peak E at 32 kev represents another radiation ; 
or alternatively, whether a radiation of about 44-9 kev energy is present 
as well, having its escape peak at 32 kev. To avoid the uncertainties 
due to escape peaks a counter containing xenon was used. With xenon 
the escape peak is 29-7 kev below the main peak, so that considerably 
less confusion is possible. 

The pulse size distribution with a xenon filled counter is shown in: 
fig. 7. The two more prominent radiations, peaks A and B are now 
seen unambiguously. Their energies are 51-8 +-0-4 key and 38-5-40-4 kev. 
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Additional radiations at 60-3-L0-4 kev (peak C) and at 45-1-L0-4 kev 
(peak D) are also present. The peak E at 69-2-+40-4 kev was investigated 
using a larger source of 78°Pu. It was shown to be the escape peak of 
a weak radiation of energy equal to 101 kev. An even weaker radiation 
of energy 115 kev was also present. The energies of the K, and K, 
radiations of uranium are 96-9 kev and 111 kev while those of plutonium 
are 101 kev and 115 kev respectively. This radiation is therefore almost 
certainly K radiation, but the measurements made were not sufficiently 
precise for deciding to which of these two elements the radiation belongs. 


Fig. 7 
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Gamma radiations from plutonium detected with a xenon filled counter. 


The plutonium samples that were used had all been carefully purified 
to remove fission products. As an additional check that the weak 
radiations observed were not due to fission products, one of the samples. 
that had been measured was put through an additional cycle of repuri- 
fication from fission products. The process was expected to reduce the 
amount of fission products present by at least a factor of ten. After this. 
additional purification process the sample again gave the same y rays 
in about the same intensities. In addition we have now observed the 
radiations with several different samples of plutonium. There is therefore 
a strong indication that the radiations are not due to fission products. 
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There remains the possibility that some of the radiations are due to 
other isotopes of plutonium or transuranic elements which are known to 
be produced at the same time as plutonium. Accordingly the radiations 
from two samples of plutonium of different isotopic constitutions were 
compared. The samples contained *3*Pu, *4°Pu and *4?Pu in small 
amounts in addition to 23%Pu. The percentage by weight of *4°Pu was 
about one half as much in sample I as in sample II. Much smaller 
percentages of 238Pu and *41Pu were present. Sample I contained about 
one third as much 2°8Pu and about one fifth as much *41Pu as sample II. 

The y-ray spectrum from sample If showed an enormously increased 
amount of the radiation of 60-3 kev energy (peak C in fig. 7) relative to 
the other radiations. This radiation was shown to be due to *44Am— 
the daughter of #41Pu (~10 years half-life) which is also present in the 
samples. The americium was removed by a chemical oxidation and 
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precipitation cycle. Both samples were treated in this way and measure- 
ments made as soon as possible after the chemical separation to avoid 
the growth of further amounts of 241Am. 

The y-ray spectra from an 8 mg source of sample I and a 6 mg source 
of sample II are shown in fig. 8. The y-rays of energy 52-0 kev and 


38-5 kev are present in both samples in nearly the same intensities and 
are therefore due to *39Pu, 
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: The 45-0-+0-3 kev radiation is about twice as intense in sample IT as 
m sample I and is therefore due to 24°Pu or possibly 2°8Pu. The intensities 
of the radiations from 78°Pu were estimated from the efficiencies of the 
counter for the two radiations. The mean track length of y rays in the 
counter was calculated for a point source on the wall of the counter. 
The mean value was 0-92 x (diameter of the counter) for y rays entering 
the counter. 

Measurements with both krypton and xenon counter fillings were used 
in deriving the energies and intensities. The final results are contained 
in table 2. 


Table 2. Radiations from 739Pu 


Estimated 
Radiation intensity per 
« particle 
Uranium L x-rays 4+1x10-? 
52-0 ++0-3 kev gamma ray lL Ome 
38:5 0-4 kev gamma ray 2x410-° 
K x-rays 2p ates 


§ 6. ORIGIN OF THE L x-RADIATION 
The L x-radiation could arise in the following ways :— 


(1) As a result of excitation produced by the passage of an alpha 
particle through the electron shell of an atom. 


(2) As a result of orbital electron capture by the nucleus. 
(3) As a result of internal conversion of a nuclear y ray in the L shell. 


We can exclude process (1) on the following grounds : since the radia- 
tion is characteristic L radiation of uranium it cannot arise from the 
passage of «-particles through the plutonium of the source. It could 
conceivably arise as a result of the excitation of the electron shell of its 
parent atom by an «-particle during its emission. This process would be 
a universal property of « emitters. Investigations so far (Rubinson 
and Bernstein 1951) suggest that the intensity of L radiation produced 
by this process is very much weaker than that observed here. Process (2) 
is energetically impossible for #°Pu since *8°Np is f--active. It is 
concluded, therefore, that the L radiation arises as a result of the internal 
conversion in the L shell of a y-ray which follows the «-particle emission 
from 239Pu. This conclusion is supported by the results of Albouy and 
Teillac (1951) using photographic plates who found soft conversion 
electrons from 28°Pu having an intensity of 0-16 conversion electrons per 
x particle. The energies corresponded to those of the L conversion 
electrons of a y-ray of energy equal to about 50 kev. The weak intensity 
of y rays observed here compared with that of the L rays or conversion 
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electrons means that there is a high degree of internal conversion (con- 
version coefficient ~ 1000). Some further experiments were performed 
to see whether the excited state of 285U in question was long lived. 
A chemical separation of the uranium daughter of plutonium was carried 
out by solvent extraction. The uranium fraction was then deposited 
on an aluminium foil and the latter placed near the window of the pro- 
portional counter. No L rays were detected and this enabled the half- 
life of the transition to be set at less than two minutes. In a further 
experiment, coincidences between « particles detected by a geiger counter 
and L rays detected by the proportional counter were measured, using a 
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coincidence resolving time of 4 microseconds. The number of coinci- 
dences observed agreed with that calculated if it was assumed that all 
the L rays observed were in coincidence with « particles. It is concluded 
that the half-life of the excited state of 235>U which gives rise to the L 
radiation is less than 4 microseconds. 


§7. MEASUREMENTS WITH 233U 
Studier (1947) has reported soft radiations from 223U, It was therefore 
decided to investigate this isotope with the proportional counter technique. 
A source containing 1:9mg of **8U was deposited on aluminium. 
L x-radiation was observed from this source by means of a krypton filled 
counter. The pulse size distribution is shown in fig. 9. The energies of 
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the peaks were measured in a similar manner to that described for 
plutonium. 


The results are given in table 3. 


Table 3 
Experimental values (kev) 12-7+0-2 16-1 0-2 19-0-0-2 
Energies of thorium L rays (kev) 12:9 15:9 19-0 
Energies of uranium L rays (kev) 13:6 16-9 20-2 
Relative intensities 1-0 1:8 0:3 


The L rays are evidently those of thorium, the daughter element in 
the decay of ?5°5U. The overall intensity was estimated as (4-41) 107 
which is closely similar to the intensity obtained with plutonium. 
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Soft y radiations were looked for with counters containing krypton 
and xenon. The pulse size distribution obtained with a xenon counter 
is shown in fig. 10. The prominent radiation at 42-8--0-3 kev is more 
intense than the radiations observed from plutonium, but is again very 
much weaker than the L radiation. The results of the measurements 
with ?°8U are given in table 4. 

By an argument very similar to that used in the discussion on plutonium 
it can be shown that the L radiation from ***U must arise from the 
internal conversion of a y ray which follows the « particle emission. 
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Only the weak y radiations listed in table 4 were found in the energy 
region up to 120 kev. A converted y radiation in this energy region must 
be responsible for the L radiation, since if the y radiation were of energy 
oreater than 120 kev K x-radiation of comparable intensity to that of 
the L radiation would have been observed. The y radiation responsible 
for the L radiation must therefore be highly internally converted. 


Table 4 
Energy of Radiations from *7U Estimated intensity per « particle 
(Thorium L x-rays) Aol x 10> 
42-8-+0-3 kev pats. 
56-1 0-4 kev 1028 


§8. CONCLUSIONS 


Proportional counters are very suitable for investigating y radiations 
in the region up to 100 kev. The presence of fine structuie in « particle 
spectra can be revealed by the presence of L radiation arising from the 
internal conversion of soft y rays. The degree of conversion is often 
very high, so that only weak unconverted y radiation is emitted. 
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SUMMARY 


Cascade theories appropriate to the various schemes of plural production 
of mesons suggested by Heitler and Janossy (and one by us) are developed. 
A set of master diffusion equations are set up which contain the various 
schemes as special cases. Their exact solution is given and detailed 
numerical results are presented for one of the schemes. These are 
presented in a form suitable for direct comparison with those obtained 
experimentally using photographic plate techniques. It is suggested, 
now that the detailed numerical results are available, that a recomparison 
be made between theory and experimental results previously obtainéd, in 
order to see whether the scheme considered does or does not agree with 
experimental observations. Results are presented in tabular form both 
for light and heavy elements and for mesons as well as protons. 


$1. INTRODUCTION 


In a series of papers (Camerini et al. 1951 a, b, 1952) under the general 
title of ‘ Nuclear Transmutations produced by Cosmic Ray Particles of 
Great Energy ’, the Bristol group published experimental results relating 
to well over 15,000 nuclear disintegrations. Their data is by far the most 
extensive to appear to date. From a detailed analysis of this data it was 
hoped to throw some light upon the fundamental processes involved in 
high energy nucleon-nucleon and nucleon-nucleus collisions. In particular 
it was hoped to give a satisfactory answer to the highly controversial 
question : Are mesons produced in high energy nucleon-nucleon collisions 
according to the multiple or plural theory of meson production ? ‘The 
Bristol group concluded that their results could be explained in terms of 
either the multiple theory with secondary generation of mesons, or the 
plural theory assuming many body collisions. They were, however, 
seriously handicapped in carrying out their analysis because of the absence, 
at that time, of detailed numerical results based on the plural theory. 


ee ee ee 
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Heitler and Janossy (1949, 1950) and Terraux (1951) have carried out 
limited calculations based on the plural theory, but these rested funda- 
mentally upon the dubious postulate that successive collisions within the 
nucleus could be treated as statistically independent events. Furthermore 
their treatment did not take the cascade character of the process fully into 
account. Detailed theoretical results based upon the various theories of 
megon production are essential for the experimentalist if he is to make 
an adequate comparison between his experimental observations and the 
predictions of theory. We shall develop in this paper the cascade theories 
appropriate to a number of schemes of meson production. proposed by 
Heitler and Janossy (1949, 1950) and give detailed numerical results for 
one of these. A further scheme, suggested by the recent experimental 
results of Lock and Yekutieli (1952) will also be developed. 

The purpose of this paper will not be to attempt to prove or disprove 
either the plural or multiple theory of meson production, but rather to 
present detailed results which will provide the basis for future comparison 
between theory and experiment. 


§2. Mmson-NUcCLEON AND Mrson-NUCLEUS INTERACTIONS 


In the original theory of Heitler and Janossy (1949) it was assumed :— 


(a) Every high energy nucleon-nucleon collision leads to the creation of 
one and only one meson. 

(6) The primary nucleon makes on the average more than one collision 
in traversing the nucleus (the effect of recoil nucleons was neglected). 

(c) The mesons produced are simply emitted without any interaction 
with the nucleons of the nucleus in which they were created. 


Later Heitler and Janossy (1950) and Terraux (1951) endeavoured to 
include the effect of recoil nucleons (hereafter referred to as scheme HJ(1)) 
and at the same time suggested two further possible theoretical schemes 
for the production of penetrating showers by nucleons and z-mesons. In 
the one scheme it was postulated that a 7-meson gives a large fraction of 
its energy to a nucleon inside a nucleus and that the nucleon then gives 
rise to a penetrating shower via the cascade process (referred to as HJ(2)). 
In the other scheme it was suggested that the 7-mesons may give rise to 
penetrating showers by a multiple process in which the primary meson 
‘ splits ° up into several mesons (referred to as HJ(3)). 

Recent experimental results of Camerini e¢ al. (1952) suggest that the 
mesons suffer only small fractional energy losses in emerging from the 
nucleus in which they were created. At the same time strong evidence 
was given which supports the view that an isolated meson suffers a large 
energy loss in traversing nuclear matter and can, when in collision with a 
nucleus, lead to the production of further mesons. The collision cross 
section for meson-nucleus collisions was shown to be geometric or very 
nearly so. These authors suggested that the apparent discrepancy 
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‘between the results might be resolved by assuming that only a few high 
energy mesons are in general created in primary collisions and that these 
mesons undergo multiplication in subsequent megon-nucleon collisions. 
Rosser and Swift (1951) have put forward evidence supporting this view. 
These authors exposed photographic plate emulsions both above and below 
30 cm of lead, at an altitude of 9000 feet, and found that for showers with 2 
‘or more ‘ shower particles’ the ratio of charged to neutral primaries is 
greater for the showers observed beneath the lead than for those above. 
They concluded that it was necessary to postulate an interaction of the 
mesons with nuclear matter whereby further mesons were created. The 
two schemes proposed by Heitler and Janossy were considered and it was 
-concluded that only the second scheme, that of meson splitting, might 
account for their results. Lock and Yekutieli (1952), in a preliminary 
study of nuclear disintegrations produced by z-mesons in the energy 
range 0:05 to 1-1 Bey, have arrived at a similar conclusion. They also 
argue that a low energy primary meson may make on the average 2 to 3 
collisions within the nucleus. 

One is tempted to deduce from the results of Lock and Yekutieli that 
very high energy mesons cascade within a nucleus in a manner entirely 
analogous to that for nucleons. We therefore now postulate that the 
mesons interact with nucleons, and that for high energies each meson- 
nucleon collision gives rise to two mesons and one nucleon, thus 


(A) meson-+nucleon—meson-+ meson- nucleon 
(B) nucleon+nucleon—nucleon-+nucleon-+ meson. 


‘The mesons and nucleons created may cascade within the nucleus and in 
dispersed matter. It is further assumed that for the high energies 
considered the cross sections for processes A and B are equal. This 
assumption is probably an oversimplification of the physical reality but 
some support is given to it by the complete symmetry of the cross section 
for process B. Ina paper by Messel (1951 a) it was shown that the best 
possible choice of F(7,, 7), the cross section for process B, was a function 
symmetrical with respect to the 7, 7, and y,,, the meson energy. This 
choice of F(y,, 72) yields average values of 71, 72 and 7,, which are equal : 
M1=N2=Nm=}- Itis reasonable to draw an analogy between A and B and 
hence to use the same cross section for both processes. It is interesting 
to note that process B has been observed in laboratory experiments, but 
process A has not, perhaps because only low 7-meson bombarding 
energies are at present available. In process B it is further assumed that 
one and only one meson is created in each nucleon-nucleon collision. 


§3. Mrxep Nucteon-MEson CascaDES 


We now give a set of master equations for the average number of nucleons 
and the average number of mesons at a depth x in homogeneous nuclear 
matter, for the various schemes of shower production proposed by Heitler 
and Janossy and by us (MPM). Let n) (H, x) dH be the average number 
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of nucleons with energies in the interval Z, dH at depth x collision units. 
in homogeneous nuclear matter, due to a primary particle (j). When 
j=1, the primary particle is a nucleon, when j=2, ameson. Leta)(H, x) 
represent the corresponding quantity for7-mesons. The master equations. 
are then given by 


{zs +1} nO”(E, 2) = (s NB x) F(E/E')dE'/E' 


ae { " )(B', 2) K(H/E')\Gw'/E', . . . . .. (Oy 
E 
(35 + aha, x)=4 | nH’, «)F( E/E’ dB’ |B" 


+3 i ” (HY, e) URE GEE’, . . . . . . (2) 
E 
where 
=J{F(n, n2)+F (ne n)}dme - - - - - (8) 


and qg, k(H/H’) and I(H/E’) are to be chosen according to the scheme 
considered. 

The solutions for the average number of nucleons N“(#) and the 
average number of z-mesons 7/)(H) with energies greater than H emitted 
from a nucleus, are obtained from (1) and (2) by making use of 
Laplace—Mellin transforms and averaging over all possible path lengths 
through the nucleus. The results are left in convenient matrix form : 


lee NG E> ae ces g—A,(s)—Us) k(s) 
aV(H)  w@(H) | * \Ag(s)—A,(s) ie 8) a(s)—A,(s) 
V[D4A2(8)] F¢—Aa(s)—Us k(s) 
X38) —Aats) LO, 8) ce welt” (%) 
where 
ite Ma-ta(s)—U(a)} * [bg +a(s)—U(s)? 
+k(s)W(0, s)—a(s){g—U(s)}}2 2. 2 . 2 2... BY 
1 8o+t0 
= 35 [oli asls oS ne 
«(s)=1—2W(0, s) rr re ac 
Lt pl—7 
Wer s=] 2M (1,2) dyy,dy2 . . (8) 


V([D,A,(s)]=2 ents ep Pahl Shad ae (9) 
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and D, is the average number of collisions suffered by a particle in making 


a diametrical passage through a nucleus of atomic weight A. 
Furthermore, 


I(s)= | ae ytelEY aye’, eee be Oh 10) 


and 
1 
k(s)= | ape ye CeB) a" pote iL) 


The values of q, k(H/H’) and 1(H/E’) are given in table 1 for the various 
schemes considered. 


Table 1 


Scheme k( E/E’) 


HJ(1)* 0 


k’ =constant 


SF(E/E’) 


In table 2 the values of I(s), k(s), A,(s) and A,(s) are tabulated for the 
three schemes. 
Table 2 


A2(8) 
0 


HJ(2) | k’/(s+1)|k’/(s+1)| 1—W0, s)—k’/2(s-+1)} 
+ {W2(0, 5) + (he’)?/4(s-+1)2}1/9] — {W2(0, 5) +(h/)3/4(8-+ 1)2} 4 


(1—W(0, s)—k’/2(s+1)} 


MPM |2W(0,s)| W(0, s) 1—W(0, 8) 1—3W(0, s) 


$4, RESULTS 


We have carried out detailed calculations for the scheme HJ(1) using 
eqns. (4) and (5). The explicit form of the cross section F(y,,72)= 
1207 2(1—7,—7,) for inelastic nucleon-nucleon collisions is taken from the 
work of Messel (1951 a), where it was shown that its use allowed one to 
account adequately for the development of the nucleon component of the 
cosmic radiation in the atmosphere. The results of the calculation are 
presented in tables 3 and 4. For the light elements carbon, nitrogen and 
oxygen an average value of D,=3-7 was used ; for the heavy elements 
silver and bromine the value D ,=6-8 was chosen. 


* Clementel (1950) has also obtained a solution for this scheme for a set of 
special boundary conditions. The solution for the nucleons has been obtained 
previously by Messel (1951 b) in his treatment of the more general problem of a 
nucleon cascade developing in the atmosphere. 
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It will be noted that in both table 3 and 4 primary nucleon energies as 
low as 1 Bev were considered. In the energy region H)>=1 to 5 Bev, the 
occurrence of elastic nucleon-nucleon collisions is important, so that our 
calculations will actually give an overestimate of the number of shower 
particles in this energy region. When using the tables, this point should 
be kept in mind. 

The symbols in the various columns denote the following (we assumed 
that one half of the nucleons were protons and two-thirds of the mesons to 
be charged) : 

P,=average number of protons with kinetic energy greater than 


500 Mev. 
P,=average number of protons with kinetic energy greater than 
800 Mev. ; 
P,=P,—P,=average number of identifiable protons among the 
shower particles. 


IT,=average number of charged 7-mesons with kinetic energy 
greater than 80 Mev. 

IT,=average number of charged 7-mesons with kinetic energies 
greater than 1100 Mev. 

IT,=IT,—I],=average number of identifiable charged 7-mesons 
among the shower particles. 

n,=l1,+P,=average number of shower particles. 

IT+ P=II,-+-P,=average number of shower particles which cannot’ 
be identified as protons or mesons because of their high 
energies. ‘ 

IT,+-P,=average number of identifiable protons and charged mesons 
among the shower particles. | 
In fig. 1 we have plotted the number of shower particles n, against the’ 
primary energy H, both for the light and heavy elements. In fig. 2 we 
give the curves for the mean square angular deviation of the shower 
protons (P,) from the direction of motion of the incident nucleon plotted 
against n,. It is likely that the angular deviation for high energy mesons 
({7,) is very similar to that for P,;. These results were taken from a paper 
in preparation by Chartres and Messel which used the theory recently 
developed by Messel and Green (1952 a, b). | 


§5. Discussion 


Even a casual glance at tables 3 and 4 reveals some rather surprising 
differences between the results for light and heavy nuclei. Though in both | 
cases the ratio of charged 7-mesons to the total number of shower particles 
ds very high, varying from 90% to 62% depending upon the energy, the 
behaviour of the various components (I7+-P)/n,, IT,/n, is very different in 
the two cases. For instance for the light elements the J-+-P particles 
constitute from 79 to 84°, and the identifiable mesons from 13 to 84. to 13% 
of the shower particles, depending upon the energy range considered. 
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This behaviour in the case of the heavy elements is completely reversed 
the IT-+-P particles now only constituting from 19 to 65%, of the shower 
particles and the identifiable charged mesons from 77 to 30%. It is 
further, interesting to note that in the case of the heavy elements there is 
steady increase in the proportion of JJ-+-P particles as the number of 
shower particles increases whereas for the light elements this proportion at 
first decreases and then increases again with increasing n,. Also for any 
given n, the proportion of JJ+-P particles is always greater in the case of 
the light elements than in the heavy. 


Fig. 1 
36 


30 
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Plot of the average number of shower particles n, against log Hp, where FH) is the 
primary energy measured in Bev. The top curve is for the heavy elements 
silver and bromine; the bottom curve for the light elements carbon, 
oxygen and nitrogen. 


The proportion of protons among the //+P particles in the case of both 
heavy and light elements is astonishingly low and in. no case does it 
increase above 40%. In the case of the light elements the saturation 
effect to be expected for high energies is readily seen ; for heavy elements 
the flattening out of the n, versus Hy curve is much slower. 

In order to make an adequate comparison between theory and experi- 
‘ment it will be necessary to apply appropriate weight factors to the results. 
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for heavy and light elements according to the proportions of these elements: 
in the photographic emulsion concerned. We have carried out this 
procedure roughly and made a comparison with the Bristol group results. 
The most surprising conclusion of this comparison has been that we were 
unable to find any great discrepancy between the theoretical results here 
predicted and the experimental results. Unfortunately, however, the 
comparison can only be made at a few points and it would be interesting 
if the experimental results could be compared again with the HJ(1) theory 
now that detailed results predicted by this theory are available. 
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Plot s the mean square angle of emission <@?> of proton shower particles: 
(P,), resulting from a nucleon-nucleus collision, against the average 
phaes " praca cia Ads The top curve is for the heavy elements 
Silver and bromine; the bottom curve for i 

( : Ur or the light elements carb 
oxygen and nitrogen. : yi 


We mention briefly one point which seems to be in disagreement with 
experimental observation, namely, the very low proportion of identifiable 
protons and high proportion of [7+ P particles. This discrepancy, if real 
would be removed in the HJ(2) and MPM models. CGilonlatigna are at 


present being carried out for the HJ(2) mod ; i 
Sie (2) model and results will be made: 
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Molybdenum, Cobalt and Tungsten at Low Temperatures 


By J. G. THomas and E. Mrenpoza | 
H. H. Wills Physical Laboratory, University of Bristol* 


[Received April 17, 1952] 


ABSTRACT 


This paper describes further results concerning the minimum of electrical 
resistance of metals at low temperatures. Some magnesium wires, known 
to contain 0-01 atomic per cent of iron, exhibited definite minima ; the 
magnetoresistance and the susceptibility were, however, independent of 
temperature. It therefore appears that though the minimum in zero field 
may well be an effect due to the presence of impurity atoms, the mechanism 
is not that of scattering of conduction electrons by unfilled electronic 
d-shells in these impurity atoms. Measurements on a dilute alloy of 
manganese in aluminium confirmed the filling of the d-shells of the 
transition metal by the multivalent solvent. 

The resistance of each of several samples of molybdenum went through 
a minimum rising to a constant value below about 0-1° kK. A supercon- 
ducting impurity (7',~0-4° K, H)~30 gauss), unidentified but probably 
a compound, could be produced by annealing. Samples of cobalt 
(hexagonal close-packed with traces of face-centred cubic) and tungsten 
showed no anomalies or superconductivity down to the region of 0-06° K. 


$1. IyrTRoDUCTION 


IN a previous paper (Mendoza and Thomas 1951) the authors described 
measurements of the electrical resistance of specimens of gold, silver and 
copper at low temperatures (i.e. from 4-2° K down to the region of 0-06° K). 
The resistance of most specimens of these metals was found to increase 
with decreasing temperature, implying the existence of a minimum in the 
resistance-temperature curve, whereas present-day theory demands a 
constant residual resistance in this region. Some of the specimens 
exhibited normal magnetoresistive effects (i.e. they obeyed Kohler’s rule) 
while others behaved abnormally. 

Similar behaviour in these and other metals has been extensively 
studied ; to the list of references in the previous paper may be added the 
following : Rorschach and Herlin (1951) demonstrated that the effect is 
still present in magnesium even when the resistance is measured without 
attaching leads to the specimen ; Serin (1951) found that the resistance of 


a specimen of gold in the helium range is independent of measuring current 
———$—$ $e See eee 
* Communicated by Dr. L. C. Jackson. 


“On the Electrical Resistance of Magnesium etc., at Low Temperatures 901 


over a range from 3 ma to 4.4. The measurements of Friedberg, Ester- 
mann and Goldman (1952) of the specific heat of magnesium seem to 
eliminate the possibility that the resistive mechanism is like that in a 
semi-conductor. The work of MacDonald and Templeton (1951) on gold 
containing added amounts of copper and nickel, and of Gerritsen and 
Linde (1951) and Gerritsen (1951) on gold, silver and copper contaminated 
with manganese, tend to concentrate attention on the role of specific 
impurity atoms in producing the effect ; the results of MacDonald (1951 b), 
whose curves for copper from different sources are different in shape, may 
imply a similar mechanism. Gerritsen and Korringa (1951) have pro- 
posed a model to explain the results of Gerritsen and Linde ; other more 
general theories have been described by MacDonald (1951 a), who relates 
the anomaly in zero field to anomalous magnetoresistance effects, and by 
Klemens (1951) in terms of the quantum theory of electrical resistance 
carried to a higher approximation. 

The present paper describes measurements on ‘ pure’ magnesium, 
molybdenum, cobalt and tungsten, and on aluminium both ‘ pure’ and 
contaminated with small amounts of manganese. The apparatus and 
methods of measurement have been described in our previous paper. 


§2. ResuLtts anp Discussion—MAGNESIUM 


All specimens used came from a batch of 0-125 mm diameter wire drawn 
by Johnson Matthey from their rod (Lab. No. 1848). A spectrographic 
analysis of the rod revealed the presence of Mn, Fe, Si, Cu, Ag, Pb, Al, 
Ca and Na, the lines being either faintly visible or very faintly visible. 
Chemical analysis of the material have the following results: Fe 0-013%, 
Mn 0-0023%, Pb 0-0013%, Cu, Al, Ag not detected. The estimated purity 
of the material was greater than 99-98%. 

Most of the experiments were done with the current and potential leads 
attached with screw contacts.* 

Resistance Measurements in Zero Magnetic Field 
A striking feature observable in the results is the great variation of the 


resistance-temperature curves from specimen to specimen cut from the 
same piece of wire. If the cause of the increase of resistance at low 


 * Of the curves taken with soft-soldered leads, we may mention specimen Mg 1 
which showed an extremely large increase of resistance with decreasing temper- 
ature. The wire was annealed for three hours at 300° ¢ in vacuo ; Tin above 
20° x; if R,=resistance at 273° K, resistance below 20° k given approximately 
by R/R,=10-2 (10-5—24/7), i.e. 0-117 at 20° x, 0-165 at 4° K, 0-265 at 15° K ; 
magnetoresistance positive at 4:-2° K, negative at 1:5°K. Later the wire was 
accidentally broken and the recoverable part was tested again; with soft- 
soldered leads there was then found only a small increase with decreasing 
temperature, with the minimum at about 4° k. These results may be due either 
to the irreproducibility of contacts soft-soldered to magnesium, or to the 
presence of a portion of the wire of high resistance and anomalous properties 
somewhere in the part which was broken off. 


> 
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temperatures is in fact the presence of certain types of impurity atoms, 
this may mean that the impurities were distributed unevenly along the 
length of the wire. 

Table 1 gives the values of 7',,;,, the temperature of the minimum, and 
Ryn» the value of the resistance at that temperature, both for our speci- 
mens and those of others whose published data quotes these values. 
There does not seem to be a universal curve connecting these two para- 
meters, at least for small values of R,,,, though there is a general trend 
linking larger 7',in with larger R,,,. In terms of the ‘contamination ; 
theory of the origin of the minima, this may be taken to imply that no 
single impurity was predominant in all specimens. (This is in contrast to 
the results for gold, where the existence of a single 7',,;,—R,,;, curve for 
specimens of widely differing origins may imply that one impurity— 
probably silver—was predominant in very many specimens.) 


Table 1 
10? X Rmin/ Ro Tin (°K) Reference 
(~100 Above 25 Present work : Mg 1) 
53 Above 20 Meissner and Voigt (1930) 
32 About 7 os a a bs 
44 16 Rorschach and Herlin (1951) 
16 5 Present work : Mg 2 
15:7 0-7 Ae : Mg 5 
8 2 & z Mg 4 
75 5 i. be Mg 3 


_ Of our specimens, Mg 2 was not annealed ; all others were annealed 
in vacuo for three hours at 300° c. 


A feature of all curves for magnesium in zero field is that the rate of rise 
of resistance increases with decreasing temperature down to the lowest 
temperatures investigated (0-06° K). 


Magnetoresistance Measurements 


The results of many magnetoresistance experiments in transverse fields 
are contained in fig. 1, a reduced Kohler diagram, using a Debye © of 
300° k. (The discrepancy with the curve given by Justi (1948) shown 
dotted, seems to be real, and not entirely due to the use of a different 
Debye @.) It is seen that over a considerable range of temperature— 


ie "ag Pk: 
0-2 K to 4:2" K—all specimens give results which are normal to within 
experimental accuracy. 


Susceptibility Measurements 


By a small modification of the Sucksmith balance, we were able to 
measure the susceptibility of a rod of the magnesium from room temper- 
ature down to 1:2". The accuracy was about 5°% at all temperatures 
Little disturbance due to the de Haas—van Alphen effect was expected x 
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the applied field was inhomogeneous. The absolute value of susceptibility 
was -+-0-25 x 10-6 per gram and within experimental error there was no 
change with temperature. Selwood (1943) quotes the same figure as the 
best value at room temperature. ; 


Discussion 

These results concerning the magnetoresistance and the susceptibility 
are of interest because the specimen was known to contain 0-06 atomic 
per cent of iron and 0-01 atomic per cent of manganese. The behaviour is 
in marked contrast to that of alloys of the noble metals containing 


Fig. 1 


4 1.25° 42° 
x 45° + 4-0° 
7 1-25° 4.0° 
@ 02° 6 1:252 
fe) 42° 


10 50 
10 -HR»/R 


comparable concentrations of iron and manganese ; these are known to 
have abnormal magnetoresistance effects (Nakhimovich 1941, Gerritsen 
and Linde 1951) and to have paramagnetic susceptibilities which vary 
strongly with temperature (Nakhimovich 1941, and see Néel 1940 and 
Constant 1945 for further references). To explain the results for the 
contaminated noble metals it is postulated that even in solution the 
transitional elements still have their d-shells incompletely filled. But it is 
also known that when transition metals are dissolved in metals of valency 
greater than unity, the d-shells may be filled up, when the atoms lose their 
transitional structure ; this has been deduced from the limits of solubility 
(Raynor 1949). The resulting alloys should then be expected to have 
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much more normal magnetic properties. Our measurements may there- 
fore be interpreted to mean that this filling of the positive holes in the 
d-shells occurs when iron and manganese are dissolved in magnesium. * 
Since however the resistance of our specimens in zero field exhibit pro- 
nounced minima, it seems that scattering of the conduction electrons by 
the unfilled d-shells of the impurity atoms—the mechanism dealt with 
by Gerritsen and Korringa—is not also responsible for the minimum effect 
in magnesium. 


r 


§3. ResuLts AND Discussion—ALUMINIUM AND ALUMINIUM-MANGANESE 


Measurements in Zero Magnetic Field 

Some ‘ super-pure ’ aluminium (99-995°%) was drawn through steel dies 
to a diameter of 0:3mm. The specimens were annealed at 300° c for 
two hours in vacuo. As the resistance was rather low our measurements 
were confined to the helium range. The residual resistance was found to 
be constant to 1 part in 10? throughout this temperature range 
(R/R = 2:89 x 10-3 for Al 1, 2-57 10-3 for Al 2). Of previous results for 
aluminium (see our previous paper for references), some have shown 
slight minima. 

In order to test the statement made in the last section—that a multi- 
valent element containing a small proportion of transition element as 
impurity may have a normal magnetoresistance—a carefully homogenized 
alloy of aluminium with 0-05°% (i.e. 0-10 atomic per cent) of manganese 
was prepared. In zero field it was found to have a resistance which was 
constant to 1 part in 10? from 4-2° K to 1-3° K (R/R,=15-5 x 10-3 for Al 3, 
15-4 10-3 for Al 4). Matthiessen’s rule was obeyed—that is, the difference 
between the specific resistance of the alloy and the pure metal was the 
same at room temperature as at helium temperature. 


Magnetoresistance Measurements 


The magnetoresistive behaviour of the aluminium and of the alloy 
with manganese was found to be normal from 4:2°K down to 1:3°x:; 
transverse fields up to 1500 gauss were used. The results are represented 
in the reduced Kohler diagram, fig. 2, which uses a Debye © of 395° x. 
(Justi’s curve, shown dotted, uses the same value.) Thus the alloy 
behaved as if the d-shells of the manganese were completely filled by 
electrons donated by the aluminium. 


§4. THe OrIGIN OF THE MinIMUM oF RESISTANCE 


We may assemble the following data : (a) the results of MacDonald and 
Templeton on the minimum shown by gold wires containing copper (and 
with them, the possibility mentioned above that silver in gold can also 


* In connection with the work on silver-manganese alloys, it may be of interest 
to remark that the curves given by Schulze (1941), for liquid air temperatures 
and above, indicate that the addition of tin, in concentrations comparable with 
that of the manganese, can remove the maxima from the resistance-temperature 
curves. Here the tin seems to act as an electron donor to the manganese, 
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produce a minimum) ; (6) the present authors’ observations on a gold wire 
with a minimum but with normal magnetoresistance ; (c) the above data 
on magnesium. These seem to show that, if the minimum is indeed an 
effect due to the presence of impurity atoms, it can be produced even by 
impurities which one has hitherto represented by models of a very simple 
kind. In dealing with the data on solutions of the transition metals in 
noble metals, quite complicated structures can plausibly be attributed to 
the impurity atoms, and these can be made to account for the experimental 
results. At present, however, it seems that comparable methods cannot 


be applied to such systems as, say, dilute solutions of copper in gold or of 
iron in magnesium. 


Fig. 2 
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§5. Resutts anpd Discusston—MoLyBDENUM 


All specimens were obtained from the B.T.H. Co. Ltd., in the form of 
polycrystalline wires, some of them 0-10 mm in diameter, others 0-025 mm. 
The principal impurities were quoted as O 0-113%, Fe 0:008% and SiO, 
0-015%, together with spectrographic traces of N iand Cr. The resistance 
of the best specimen was lower than any previously recorded except that 
of Justi and Foroud (1939) ; but as their specimen was a single crystal and 
ours a polycrystalline wire, the chemical impurity contents of the two 
specimens were probably about the same. A . 

Leads were attached by silver soldering using a nichrome soldering bit. 
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Resistance Measurements 

The results for several specimens are shown in fig. 3. Details of the 
specimens and of the measuring conditions are as follows: Specimen 
Mo 3: diameter 0-025 mm, no heat treatment, measuring current 0-06 
and 0:1 ma. Specimen Mo 4: next to Mo 3 on the reel of wire ; with no 
heat treatment gave a curve close to that for Mo 3. Later annealed for 
two hours in a vacuum better than 10-* mm Hg, most of the wire being at 
800° c, with some ‘ hot spots’ at 1200°c ; measuring current 0-2 ma for 
the lower branch of the curve. The fall of resistance was later realized to 
have been due to the presence of small amounts of a superconducting 
impurity ; the curve was taken again in the presence of a magnetic field 
of 30 gauss using 0-2 mA measuring current ; the upper branch was then 
obtained. The superconducting impurity will be dealt with later. 
Specimen Mo 2: diameter 0-1 mm, no heat treatment, measuring current 
0-3 ma with check points in the helium range using 3maA. During part 
of the work in the helium range, measuring conditions were not good. 
Specimen Mo 1: diameter 0-1 mm, no heat treatment, measuring current 
0-3maA. Specimen Mo 5: diameter 0-1 mm, annealed at a uniform 
temperature of 1 500° c for two hours in a vacuum better than 10-> mm Hg; 
measuring current 0-3 ma for the lower branch—which went to zero 
resistance—the upper branch being taken partly in zero magnetic field 
using 3mA measuring current, partly in fields up to 30 gauss using 
0-2 mA measuring current. 

A notable feature of the curves for the normal metal is the constancy of 
resistance at the lowest temperatures. This is thought not to be a spurious 
effect for the following reasons: The flattening was observed under 
conditions where, with other metals, the resistance was found still to rise 
with decreasing temperature. Further, the curves were reproduced within 
close limits using a range of measuring currents, and the flattening is 
therefore not thought to be an effect due to poor thermal contact between 
the salt and the wire. Lastly, whereas most of the points above 0-05° K 
were obtained by demagnetizing to the region of 0-05° K and allowing the 
specimen to warm up slowly, others were obtained by direct demagneti- 
zation from fields less than the maximum available ; both sets of points 
gave consistent results and showed that no appreciable persistent temper- 
ature inhomogeneities (which might have produced errors) were present in 
the salt. 

A further feature is that over a considerable range of temperature the 
resistance varies exponentially with decreasing temperature. None of the 
curves given by Gerritsen and Linde exhibits either the flattening over the 
2: 1 temperature range or the exponential variation. 

It may be noted that Meissner and V oigt (1930) found an increase of 
resistance at low temperatures, Steiner and Fiinfer (1936) found none. 
The Superconducting Component 

The curves for Mo 4 and Mo 5 showed that a superconducting component 
was present. ‘The superconductivity could be destroyed either by applying 


Resistance of Magnesium etc., at Low Temperatures 907 


108 


106 


104 
94 
EOS 


982 


90 


88 


x10 


Temperature (deg. absolute) 


SER. 7, VOL. 43, NO. 343.-—AUGUST 1952 3a 


908 J. G. Thomas and E. Mendoza on the Electrical 


a transverse magnetic field or by passing a current of a few milliamps 
through the wire. At 0-3°k the values were 18 gauss to restore half 
resistance and 26 gauss to restore full resistance with 0-3 ma current ; 
with no external field 3 ma restored full resistance. Comparison of these 
two effects suggested that the superconductor was in the form of fine 
filaments and this was in accord with the great breadth of the transition in 
zero field. Presumably the filaments were continuous in Mo 5 but not 
continuous in Mo 4. 

The superconductor was only present in those wires which had been 
annealed. It has not been possible so far to identify the impurity, but 
x-ray examination of Mo 5 and a length of unannealed wire did confirm 
the presence of a second component in Mo 5. Since the position of 
molybdenum in the periodic table suggests that it might be a super- 
conductor, it is relevant to enquire whether the annealing above 1 200° c 
could have produced another crystalline modification. The x-ray work of 
Edwards, Speiser and Johnston (1951) eliminates this possibility, however. 
It appears therefore that the superconductor must be a compound of some 
sort ; for previous work on the superconductivity of molybdenum see 
Shoenberg (1940). 

Magnetoresistance Hxperiments 

Specimen Mo 3 was abnormal in its magnetoresistive behaviour ; at 
4:2° k the effect was small and positive, at 1-2° kK the value of 4R/R went 
through a maximum with increasing field, then fell to zero and became 
negative at high fields. All other specimens followed Kohler’s rule in 
transverse fields, as may be seen from the reduced Kohler diagram, 
fig. 4, using 9=450° k. The discrepancy with the curve given by Justi 
(1948) may be due to the use of a different value of 0. 


: §6. RESULTS FOR COBALT AND TUNGSTEN 
Cobalt 


Cobalt wire 99-95°/, pure was drawn for us by New Metals and Chemicals 
Ltd., to a diameter of 0-13 mm. For an unannealed sample, which was 
very springy, the value of R/Ry was 0-162 at 4:2°K, 0-160 at 1-20° x, 
remaining at that value down to 0-10° k. For a specimen annealed for 
three hours at 1000° c in vacuo, the value was determined as 0-0302 from 
4-2° k down to 1-36° k, and 0-0303 from 1-20° K down to 0:06° K. It may 
be noted that both of the specimens of Meissner and Voigt (1930) showed 
definite minima. 

An X-ray examination of the untreated wire showed a hexagonal close- 
packed structure ; the heat treated wire contained also a small amount 
of face-centred cubic structure, so that the absence of superconductivity 
for both these modifications has been shown down to 0-06° K. 

The magnetoresistance effect was the same at 4:2°K as at 1-2° K, 
showing a marked saturation ; but since in these experiments the wire was 
wound on a vertical cylinder with the field applied horizontally, the field 
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was partly transverse and partly longitudinal, and there is little quanti- 
tative value in the results. 


Tungsten 


For a single specimen of tungsten wire (supplied by the B.T.H. Co. Ltd., 
impurities SiO, and Al,O, each 0-025% ; prepared by sintering) 0-025 mm 
diameter, the residual resistance was large (R/R,—0-10) and independent 
of temperature from 4:2° kK down to 0-10°K. Other recent work on the 
same metal has given similar results (van den Berg 1948). 

The magnetoresistance effect was small but normal. 


Fig. 4 


Motrenl- 25. 
Mo2 © 1:25° 
Mo4. x 1-2°4.2° 
MoS 4 4-2° 


§7. CONCLUSIONS 


In addition to the conclusions presented in § 4 it has also been found that 
several specimens of molybdenum present further examples of resistance- 
minima with no correlated magnetic anomalies, within experimental error. 
The resistance in zero field becomes constant again below an attainable 
low temperature, whereas the resistance of all the other metals studied 
continues to rise at the lowest temperatures. 
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Ice Crystals of Spiral Form Grown from the Vapour 


By B. J. Mason 
Department of Meteorology, Imperial College, London 
and 


P. G. Owston 
I.C.1. Research Fellow, Department of Chemistry, University College, London* 


[Received May 13, 1952} 


THIs note describes some ice crystals of unusual spiral form found 
growing on the cooling pipes of a cold store containing whale-meat. They 
had been growing for about three and a half years at a constant air 
temperature of —10-+0-1° c. 

The crystals have much in common with those described as ‘ crevasse 
hoar’ by Seligman (1936), and occur in remarkably light skeletal masses 
made up of many incomplete hexagonal prisms with ‘ hopper’ develop- 
ment of both prism faces (rectangular) and basal faces (hexagonal) as 
shown in figs. 1 and 2 (Plate LX). These seem to have been formed 
by the addition of rod-like crystals of thickness } to 1 mm to the edges 
of a small central hexagonal prism. 

None of the hexagonal basal faces, however, was complete, but showed 
the spiral formation called ‘scrolls’ by Seligman, the outer turns of 
which exceeded two inches in diameter for the largest crystals. A number 
of single spirals are shown in fig. 2, and a particularly beautiful and 
symmetrical double scroll in fig. 3 (Plate LX). 

The directions of the crystal axes in different crystals varied and bore 
no obvious relation to thermal gradients or air currents. X-ray 
examination showed the rod-like crystals to be mono-crystalline and 
elongated parallel to one of the crystallographic a axes. The hopper 
faces are therefore (1010) and (0001). 

No crystal was observed with both basal faces developed. Crystals of 
various spatial orientations were detached and tested for piezo-electric 
properties by the Giebe-Scheibe method; no indication of these 
properties was found either before or after light crushing of the crystals 
to break up any possible twins. Pyro-electric properties were also not 
found though the apparatus used could detect charges only one-twentieth 
as strong as those developed by the very weakly polar s-triphenyl 
benzene. Development at one end only of the ¢ axis is therefore not proof 
that ice possesses a polar structure. 


* Communicated by the Authors. 
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Though the crystals are of spiral form and sometimes exhibit pairs of 
spirals of opposite hand (figs. 2 and 3), this is not an example of crystal 
growth by spiral dislocations, for the crystals are made up of successive 
two-dimensional spiral layers, not of one continuous three-dimensional 
spiral. A possible explanation of the mode of formation is that at an 
early stage in growth some obstruction—perhaps the support on which 
it grew—prevented the crystal from completing a hexagonal edge and so 
initiated the spiralling; this would be consistent with the occasional 
occurrence of symmetrical double spirals. 


One of us (B. J. M.) wishes to thank the Superintendent of the Low 
Temperature Research Station, Cambridge, for permission to work in 
the cold chambers and both authors are indebted to the staffs of their 
departments for photographic and other assistance. 
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The Extra Electrical Resistance due to Cold Work and Neutron Irradiation 
of Platinum 


By R. A. DugpALE 
Atomic Energy Research Establishment Harwell, Didcot, Berks.* 


| Received June 5, 1952} 


In a recent paper, Seitz (1952) has put forward the hypothesis that part of 
the extra resistance brought about by the cold working of metals is due 
to the creation of vacant lattice sites by moving dislocations. In this 
connection it seems worth while to report the results of some recent 
experiments on the annealing of the extra resistance produced in platinum 
wire by cold work and by neutron bombardment. 

Neutron bombardment in a pile may be expected to form localized 
groups of vacant lattice sites and interstitial atoms, since, in general, a 
fast neutron on collision with one of the atoms of the lattice will give it 
sufficient kinetic energy to displace several other nearby atoms (Seitz 1949). 
On subsequent heating one may expect the vacant lattice sites or the 
interstitial atoms, or both, to acquire some mobility, leading to one or 


more of the following processes :— 
(1) “ Recombination ’ of interstitial atoms and vacancies. 
(2) A randomizing of their spatial distribution. 

(3) Condensation of vacancies. 


4) Disappearance of vacancies and interstitials at grain boundaries. 
eS ee ee 
* Communicated by the Author. 
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These processes would lead to a reduction in the excess resistance which 
had been produced by the neutron bombardment. If a part of the excess 
resistance caused by cold work is due to lattice vacancies, one might 
expect some similarity in its behaviour on annealing. 

The decay of the extra resistance of two irradiated and one cold worked 
resistor of platinum (each of about 1 2) has been followed part of the way 
when they were heated first to 70° c for 640 hours, and subsequently to 
90° c for 336 hours. The resistance was measured at intervals to an 
accuracy of | part in 10°, after cooling to the triple point of water. The 
two irradiated resistors had increased in resistance by 0-2°%% during 
a three months irradiation at 50° c in the Harwell pile B.E.P.0. The 
extra resistance of the cold worked resistor amounted initially to about 
15%. During the two annealing periods about 13% of the extra 
resistance of the cold worked resistor and one irradiated resistor, and about 
17% of the extra resistance of the other irradiated resistor disappeared. 

All the resistance decay curves could be fitted to an empirical formula 
of the type 

ror 1 10g eae) ke at (L) 


where is the extra resistance, and A, B, and y are constants at constant 
temperature. The fit was good everywhere except in the early stages of 
the annealing at 70° c, where a small deviation was noticeable. It should 
be emphasized that this empirical relation between 7 and ¢ cannot hold 
indefinitely, and no physical significance can be attached to the constants. 

Although it is not known which of the four processes postulated above 
is the important one at these temperatures, an activation energy can be 
estimated if we assume 


dr 
di =f (r)e— OM pubed Same cars 4! « (2) 
The gradient d7/dt can be determined for any value of 7 from the empirical 


formula (1). If this is done at the point where the temperature was 
changed, U can be determined from the equation 


ono ()./@QaI/GR-# -- ° 


This has been done for all three resistors, and the following table shows the 
values of U obtained, together with the constants A, B, and y. 


Resistor A (uQ) B (pQ/hr.) y (hrs.) U (ev) 


70°c 90°c | 70°c 90°c | 70°c 90°C 


. Irradiated 9748 3592 182 570 67 250 
. Irradiated 2582 2913 146 332 56 178 
3. Cold worked 15926 15317 928 1050 36 85 
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The agreement in the value of U for the three resistors is noteworthy. 
This, together with the general similarity in shape of the annealing curves, 
suggests strongly that the same processes are responsible for the decrease 
of resistance in the two cases. This supports Seitz’s speculation that the 
increased resistance due to cold work is at least in part due to the creation 
of lattice vacancies (and perhaps interstitial atoms). 


The author’s thanks are due to Professor N. F. Mott for drawing his 


attention to the paper by Seitz, and to Sir John Cockcroft, Director of the 
A.E.R.E., for permission to publish this note. 
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Cloud Chamber Photographs of the Cosmic Radiation. By G. D. RocHESTER 
and J.G. Witson. [Pp.128.] (London: Pergamon Press Ltd.) £3 10s. Od. 


Aw atlas of cloud chamber photographs of various phenomena associated with 
Cosmic Rays is very welcome. The photographs are collected from laboratories 
throughout the world, and include a large number that have not been published 
hitherto. With each photograph there are captions which explain most ade- 
quately the main features of the event, and also stress any doubt there may 
be to the correctness of the interpretation put forward. The first section is 
devoted to technical features of operation. In the remaining sections the 
photographs show, in an extremely vivid manner, a wide variety of cosmic ray 
phenomena ranging from electron cascades to the recently discovered V 
particles. Such a book is very dependent upon the quality of production ; in 
this case it is excellent. 

The book besides being of great assistance to those who work with cloud 
chambers should in addition have a very wide appeal, since it depicts many 
fundamental phenomena clearly. Pits 


Low Temperature Physics. By F. E. Stwon, N. Kurt, J. F. Auuen and 
K. Menpetssoun. [Pp. 132.] (London: Pergamon Press Ltd.) Price 21s. 


Tne four lectures which are reported in this volume are by Professor F. E. 
Simon (A General Survey), Dr. N. Kurti (The Temperature Range below 1° x), 
Professor J. F. Allen (Liquid Helium) and Dr. K. Mendelssohn (Superconduc- 
tivity). The course was given at the Royal Institution in March 1950, but the 
version for the press contains references to work done as recently as a few 
months ago. 

This book provides an excellent account of the present position in a rapidly 
changing subject, reviewing its topics against the background of physics as a 
whole and putting the successes and the difficulties into perspective. It should 
interest not only specialists in low temperatures, but also workers in the many 


other branches of physics to which low temperature techniques are being 
increasingly applied. EK. M. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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The area of interest lies around the The lower arrow marks the layer edge 
two dislocations near the point under discussion, while the hori- 
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The hook-shaped layer edge lies The vight-hand pair of horizontal and 
near the point of intersection vertical arrows mark the hook-shaped 
of the arrows. X 625 layer edge. The left-hand pair of 
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If the three layer edges which are arrowed The area around the point of intersec 
are traced upwards it will be found of the arrows is that under discus: 
that all three finish with a hook-shape. x 
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Fig. 16 
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The three layers discussed are marked by arrows. The two twin dislocation loops lie arow 
If the two upper layer edges are traced down points of intersection of the two v 
it will be seen that they end on dislocations. arrows with the horizontal arrow. 
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The two layers with short straight edges close to the dislocations lie near the 
point of intersection of the vertical arrow with the lower horizontal arrow. 
The higher horizontal arrow marks a layer edge at the termination of which 
the dislocation has moved upwards. x 625 


Fig. 18 
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The three dislocations which have moved lie at the termination of the three 
layer edges which are arrowed. X 625 
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The sections of the layers discussed lie near 


the point of intersection of the arrows. 625 
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The vertical arrows denote the line along which the dislocation has moved, and along whick 
linear kinking can be seen. The horizontal arrows denote the limits of the area in which 
location has moved. To help in the identification of the area it is shown isolated on ther 


A limited slip-zone, formed just before, or after, growth ceased, in which it is 
possible to trace the kinking of the outer layers. x 625 


“L. J. GRIFFIN. Phil. Mag. Ser. 7, Vol. 43, Pl. LVII. 


UA Ad | 
— 


KEE EUG CU ee Ee. Vy 
See Vek AAA 
< Ayah AA 
ye \ AXA Waa’ 


Ly 
* *y 


KAA 


x 625 


Limited slip-zones, formed just before, or after, growth ceased in which it is 
possible to trace the kinking of the outer layers. x 625 
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A small-scale mosaic of a complex * buckled’ strip which has probably been produced 
by the simultaneous action of a number of different dislocation sources. x 250 


[04108 9Tqnop [eVoLAeULUIAS Vv JO SMOTA OAT, 
(‘SSOde SOYOUL OMY ULYY oL0UL ore speatds ysoSae] Jo suINy 109nQ) 
“doquuvyo plod eyy fo sedid surjooos oy} uo ns ws Surmoss speyshiQ *Z% “SLT 
‘soovy wistid pue [eseq yjoq Jo yueutdoyeaop 


2 


_leddoy, surmoys susiud yeuosexey a 


x 
= 
a 
= 
Oa 
Ke) 
> 
~ 
ce 
o 
n 
bo 
i) 
E 
= 
a 


G. OWSTON 


MASON & P. 


Bar 


